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Abstract 

In this paper a stochastic reaction diffusion system is considered, which models the spread of a finite popu- 
lation reacting with a non-renewable resource in the presence of individual based noise. A two-parameter phase 
diagram is established to describe the large time evolution, distinguishing between certain death or possible life 
of the population. 
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1 Introduction 

1.1 Statement of results 

In this paper we study the stochastic reaction diffusion system on R d , for d < 3, 

dtu = Au + ftuv — 71* + \fu W, 



(1-1) 

OtV = —uv. 

' One interpretation for this equation is that the solution u > describes the distribution of a population on R d 
r£ , at times t > and that v > is the density of a nutrient gradually used up by the population. The noise W 
C$ ' is a space-time white noise, and the multiplicative factor y/u models 'individual based' noise. This noise arises 
in situations where each individual in the population contributes independently to the noise, so that the variance 
of the noise is proportional to the population density. In dimensions d < 3, linear scaling reduces the number of 
possible parameters to two. We have chosen to take the diffusion coefficient, the rate at which the nutrient is used, 
and the noise coefficient all to be one, leaving two parameters ft, 7 > 0, which we think of as a reaction rate and a 
death rate. 



The main results concern the long time behaviour of solutions. Without noise there are travelling wave 
solutions and it is expected that sufficient noise can destroy these waves. For our description of the long time 
behaviour, we take an initial condition uo — fj, in A4(R d ), the space of finite measures in R d , and an initial nutrient 
level vq(x) taking a constant value which, by scaling again, we have taken to be 1. We say that certain death occurs 
for the parameters values ft, 7 if 

P [ut = for large t] = 1 for any initial condition fi E A4(R d ) and vq = 1. 
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We say that possible life occurs for the parameter values (3, 7 if 

P [u t 7^ for all t] > for any non-zero initial condition jj, € A4(R d ) and Vq = 1. 



The main results of the paper describe a phase diagram for the parameter values /?, 7 at which certain death or 
possible life occur. 

Theorem 1. Consider solutions to (1.1)) with initial condition vq = 1. Then for any values of the parameters 
(3, 7 > either possible life or certain death occurs. 

• When d—3, there exists a non- decreasing function (3 — > *5f((3) € (0, 0) so that when < 7 < ^{{3) possible life 
occurs and when 7 > ^((3) certain death occurs. Moreover 



• When d—2 there exists a critical curve [3 — > ^(f3) £ [0, /?) as in dimension d = 3 (but allowing the possibility 
that *(/3) = /or small (3). 

• When d—1 certain death occurs for any 7, (3 > 0. 
Remarks 

1. The estimates in the paper imply certain bounds on the critical curve. In section |6~31 we show that \&(/3)//3 — > 1 
as /3 — > 00. This asymptotic is not sharp, and the correct large (3 asymptotic is, we believe, (3 — ^{(3) ~ (3~ 2 /( 6 ~ d \ 
The remarks in section 16.31 indicate that a certain improved spatial Markov property is required to establish these 
asymptotics, and the details are left to a subsequent paper. Many other questions remain, such as regularity of 
W and the behaviour on the critical curve. Most importantly perhaps, we do not know the small (3 behaviour in 
d = 2, though we conjecture (see section 1773")) that there is certain death for 7 = and small enough (3. 

2. The possible life behaviour is not an example of global coexistence of two species, for locally either the population 
ut becomes extinct or the nutrient levels converge to zero (see Lemma [9]). Rather one expects the process to live 
on a moving front travelling through space, and behind this moving front the nutrient is used up and the process 
dies away. It is therefore more analogous to the weak survival results for the contact process on a homogeneous 
tree (see [5] section 1.4). 

3. A key difficulty in dealing with many reaction diffusion systems is the lack of pathwise comparison results 
(which were crucial to the arguments for the scalar equation studied in [3]). For example if two solutions u,u 
satisfy u < u one should not expect that u t < u t for t > 0. To replace these arguments, the key new method in 
this paper is the use of comparison results for the total occupation measure J °° Utdt and total exit measure u^P of 
solutions on a domain D. Some intuition behind this is given in section fl. 21 

4. Using scaling (as in Lemma |3]) one can investigate the dependence on other parameters. For example the 
equation 



can be reduced to the standard form (II. lj) . by a linear change of variables, with the parameters becoming 7 = 
7oo-~ 2 /( 4-d " ) and [3 — a~ 2 /( 4 ~ d ). The main theorem shows that possible life occurs, for example when d = 3, if 
< 70 < &(o~) for some $ : (0, 00) — > (0, 1) and satisfying lim^o ^(° P ) = 1 an d linv^oo ^(V) = at certain rates. 
One reason for choosing /3, 7 as our parameters is that we have natural monotonicity in these parameters which is 
unclear for other choices - for example we do not know if <!> is decreasing. 

5. One can investigate, via phase plane methods, the possibility of travelling waves for the corresponding deter- 
ministic system dtu = Au + uv — ju, dtv — —uv with vq = 1. (With no noise term and (3 > we can remove 



< liminf /3" 2 ^(/3) < limsup /3~ 2 *(/3) < 00, 
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one more parameter by linear scaling, leaving only 7 > 0.) We require 7 < 1, else the death term exceeds the 
reaction term resulting in death. Fix k £ R d with |fc| = 1. Looking for a travelling wave solution of the form 
Ut(x) = U((x ■ k) — ct), v t (x) — V((x ■ k) — ct) with c > 0, we set W = U' to obtain the first order system 

U' = W, V' = -UV, W' = -yU-cW-UV. 
c 

The steady states are (U, V, W) = (0, a, 0) for any a. By scaling we assume the untouched nutrient level is 1 and 
seek a path in phase space connecting (0,0,0) to (0,1,0), so that u is a travelling 'hump', behind the wave the 
nutrient is used up and ahead it is untouched. Linearizing around (0, 0, 0) we find the eigenvalues satisfy A = 
or A 2 + cA — 7 = 0. So for any 7 > we have real eigenvalues Ai < = A2 < A3. Linearizing around (0, 1,0) 
we obtain A = or A 2 + cA + (1 — 7) = 0. For 7 < 1 we get Ai < A 2 < = A3 when c 2 > 4(1 — 7). Complex 
roots, which imply the impossibility of a (non-negative) travelling wave, occur when c 2 < 4(1 — 7). The equations 
are three dimensional and degenerate and we have not done a rigorous analysis. Based on the above local picture, 
the following seems the reasonable best guess: there is a unique path connecting (0,0,0) to (0,1,0) which stays 
in the good region U, V > precisely if c 2 > 4(1 — 7) and hence a travelling wave for all 7 < 1 and with speeds 
greater than or equal to 2-^/1 — 7. The main result supports the conjecture that the deterministic travelling waves 
are stable to small enough individual based noise noise in dimensions d = 2, 3. 



1.2 A description of key techniques and the layout of the paper 

In this section we explain some of the intuition and tools used in the proofs, emphasising the new ideas. 

When (3 = 0, solutions u to (jl.ip are Dawson- Watanabe processes with underlying Brownian motion and mass 
annihilation at rate 7. In dimensions d > 2, Dawson- Watanabe processes are singular measures (u t (dx) : t > 0). 
However, as shown by Sugitani [13], when d < 3 the occupation measures f*u r dr are absolutely continuous with 
repsect to Lebesgue measure, with densities u(s, t, x). This will remain true for for solutions to (jl.ip and hence 
the the formal solution Vt(x) = vq(x) exp(— u(0, i, 2)) will make sense as a function. Using this idea, we construct 
solutions to (jl.ip in d < 3 by a change of measure starting from a Dawson- Watanabe process; a reverse argument 
yields uniqueness in law of solutions. 

The results on possible life use a construction of a supercritical oriented percolation process, coupled to a 
solution of (jl.ip . so that if the oriented percolation survives so too does the solution. The use of an embedded 
oriented percolation to show survival was first developed for interacting particle systems (see Durrett [1] for many 
examples). Similar methods have been used before by the authors to establish a one-parameter phase transition 
for a noisy version of the KPP equation [5]. The results on certain death in dimensions d > 2 are simpler, and the 
background idea is a comparison with a subcritical branching process (for which simple moment estimates show 
death) . 

When possible life occurs, one expects that the process lives on a front propagating with asymptotically linear 
speed. Noisy oscillations around this linear speed mean that it is not possible to predict exactly where in space, 
at a given time t, to look for the living part of the solution (that is where it is in the process of exploiting the 
nutrient). This leads to the wish to apply the percolation comparison in space at random times. To overcome 
this technical problem, we freeze the process at random stopping times, by means of exit measures. The exit 
measures can be constructed by solving (jl.ip on a domain D with Dirichlet boundary conditions and measuring 
the total flux of mass out of the boundary of the domain. Suitable comparison results do hold for this model, not 
for u t {dx) but for the total exit measures u®£{dx). For example if uo(dx) and uo(dx) are supported inside D and 
u Q (dx) < uo(dx), then one can construct solutions for which u^(dx) < u^'(dx). It turns out that the total exit 
measures are stochastically monotone both in the initial conditions and in the parameters /?, 7. A spatial Markov 
property, analogous to that found for Dawson- Watanabe processes by Dynkin in 5 , shows that exit measures can 
be constructed iteratively on larger and larger domains, and is used to set up the percolation comparison. It is 
this spatial Markov property especially, that means the use of the special form of noise ^JuW is needed for our 
methods. 
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Intuition as to why comparison results hold for the total exit measures is easiest for particle approximations 
to (jl.lj) . where population particles may react with nutrient particles leading to death of a nutrient particle and 
the creation of extra population particles. The total exit measures do not depend on the time at which particles 
hit the exit of the domain D. It does not affect the final level of the nutrient Voo{x) if we change the order that 
particles pass through position x. For example, we may think of the initial mass uo{dx) as being composed of blue 
particles, identical to those from u (dx) and some extra red particles represented by uo(dx) — u (dx). We freeze 
the extra red particles until all the blue particles have died or are frozen on 3D. Then we may run the red particles 
through whatever is left of the nutrient, but this can only lead to a larger exit measure. The reader may wish to 
draw a diagram illustrating the possible interactions for the simplest non-trivial case of two population particles 
and a single nutrient particle - it was such a diagram that initiated this work. 

Layout. In section [2] we state the existence, uniqueness and comparison results for occupation and exit 
measures. Section [3] contains the change of measure arguments needed for the proof of existence and uniqueness, 
and the proof of the existence of a critical curve ^(/J). Approximations to solutions, using a discretized nutrient, 
are given in section [H together with tightness and convergence results. Comparison results for the approximations 
are established in section [5] together with the weak convergence arguments allowing us to pass to the continuous 
limit. However, the arguments for life (in section[6|) and death (in section[7J can be read after the statement of the 
comparison results in section [2l 



2 Statements of existence, uniqueness and comparison results 
2.1 Definition of solutions 

We need to consider solutions to (jl.ip on domains D. Throughout the paper our domains are taken to be either (i) 
connected open sets with smooth boundary; or (ii) open boxes, that is of the form D = Hi( a '' f° r ~°° — a » < ~ 
hi < oo. This restriction allows us to quote estimates on the Green's functions for the domain. Moreover, boxes will 
be all that we need for the blocking arguments in sections [6] and [Jj We will look for solutions with zero Dirichlet 
boundary conditions, that is where the population is killed on reaching the boundary. This will be encoded by the 
test functions for the martingale problem below. We let C£(D) be the space of fc-times continuously differentiable 
functions <fi : D — > R for which <fi and its derivatives can be extended to be continuous bounded functions on the 
closure D. Let Cq(D) be the subset of those </> <E C%(D) for which = on 3D. Note that solutions on D = R d 
are a special case of the definition below, and in this case Cq(D) reduces to the space of bounded functions with 
two bounded continuous derivatives. 

Notation. We write A4(D) for the space of finite measures on D with the topology of weak convergence. For 
fj, £ M and (suitable) measurable functions /, g : D — > R we will write //(/) for J D f(x)fj.(dx) and (/, g) for 
Id f(x)g{x)dx. 

We look for solutions (u,v) = (ut(dx),vt(x) : t > 0, x € D) to defined on some filtered probability 

space (fi, {Tt)iT,P). The solutions will be in the sense of martingale problems. The adapted process u will have 
continuous paths in Ai(D). Let V be the [Tt] predictable sets in £1 x [0, oo) and B{D) be the Borel subsets of 
D. The process v will be predictable and non-increasing, that is the map (ui,t,x) — > v t (x)(u>) will by V x B(D) 
measurable, and such that the paths t — > v t (x) £ [0, 1] are non- increasing for all x, almost surely. We say that such 
a pair (u, v) is a solution to (|l.ip on D if for <fi e Cq(D) 




P a.s. 



(2.1) 



(2.2) 
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where mt(4>) is a continuous (Ft) local martingale, with mo((f>) — OP a.s., and with quadratic variation 

[m(4>)] t = [ u s ((/) 2 )ds, P a.s. (2.3) 
Jo 

If, in addition, P[uo — /i, vo = f] = 1, for some measurable / : D — > [0, 1] and /i € .M(Z?), we say the solution has 
initial condition (/i, /). 

We will wish to keep track of the flux of mass that exits to the boundary dD by time t, which creates the so 
called exit measure uf D (dx) on dD. For a deterministic equation, this flux would be given by a surface integral 
Jo IdD ^ U ' NdS ds for an outward unit normal N. A simple way to treat this process in the noisy setting, where 
these derivatives do not exist, is to extend the martingale problem formulation to test functions (f> £ Cj^(D). The 
exit measure will take values in M(dD), the space of finite measures on dD with the topology of weak convergence. 
We look for an adapted M.(dD) valued process t — > uf D , with non-decreasing continuous paths, so that Uq D = 
P a.s. and for £ Cj;(D) 

u t ((f)) = uo (<t>) + / u s (A<f> + /3v s (f> -j</})ds-Ut D ((/))+ m t (4>), P a.s. (2.4) 
Jo 

where mt(4>) is a continuous (Ft) local martingale, with mo(4>) = P a.s., and with quadratic variation as in (|2.3p . 

The following existence and uniqueness result is shown in section 13.21 using a change of measure argument 
starting with a Dawson- Watanabe process on D. 

Notation. We let C([0,oo), E) be the space of continuous functions with values in a metric space E, with the 
topology of uniform convergence on compacts. Both Ai(D) and A4(dD) are metrizable as Polish spaces. Let H,£> 
be the space C([0, oo), Ai(D) x M(dD)). We write (U t , U^ D ) for the canonical random variables on fi^, U for the 
Borel subsets of H,£> and (U t ) for the canonical filtration. Note that in the case D = R d where dD = the second 
component of the path space becomes trivial. 

Notation. Let B(D, [0, 1]) be the space of Borel measurable / : D — > [0, 1]. We give B(D, [0, 1]) the sigma field 
generated by (f,I(A)) for all bounded Borel AC D. 

Theorem 2. Suppose d < 3 and D is a domain as described above. Fix 0, 7 > 0. 

(i) For any /j, £ Ai(D) and f G B(D, [0, 1]), there exists a solution (u,v) to on D with initial conditions 

(ii) For any solution, the law of (u t : t > 0) on C([0, 00), M(D)) is uniquely determined by the law of (uq,vq). 

(Hi) For any solution and < s < t, the occupation measure f*u r dr is absolutely continuous with respect to 
Lebesgue measure, and there is a continuous version of the density (s,t,x) — > u(s,t,x) on < s < t, x 6 D , 
almost surely. Moreover, setting u(0,t,x) — lim s |o u(s, t, x) we have, 

v t (x) = v Q (x)e~ u ^' t ' x) for all t > 0, for almost all x £ D, P a.s. 

(iv) For any solution (u, v) there exists an adapted exit measure process u — (uf D :t>0), with non- decreasing 
continuous paths in M.(dD), satisfying the extended martingale problem \2.J$ . Moreover, for all t > 0, there 
exists a measurable map R t : C([0, t],M(D)) -> M(dD) so that uf D = R t ((u s : s < t)) almost surely. 

(v) Let Q^'f' 1 be the law of (u, u dD ) on £lr> for a solution with initial condition (fi, /). Then (/!, /) — * Q^^' 1 (B) 

is measurable for any Borel B C n D . Moreover, the family of laws (Q^j' 1 ■ V € M(D),f £ B(D, [0,1])) 
satisfies the following strong Markov property: for any solution (u,v) defined on (f2, (Ft), T, P), and any 
finite (Ft) stopping time t, 

E [h(u T+ .)\F T ] = Qf;j7 exp( _„ (0 ,.)) UU)] P a.s, 
for all bounded measurable h : C([0,oo), M.(D)) — > R. 
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We note one consequence of these results: using part (iii), whenever (u, v) is a solution we may replace v by 
the process Vo(x) exp(— w(0, t, x)) and it will still be a solution. 

The following scaling lemma will be very useful in allowing us to transfer large or small parameters to terms 
in the equations that are convenient. We hnd our intuition stronger when working in regimes with small or large 
values of parameters, rather than small or large space and time scales. This lemma allows us to transfer between 
these two regimes. 

Lemma 3. Suppose that (u, v) is a solution to il.ll) on D. For a, b, c, e > define 

u t {A) = ^-tum(cA), v t {x) = evbt(cx) 
c a 

where cA = {cx : x £ A}. Then (u,v) is a solution on c~ 1 D, in the sense described above, to the equation 

d t u = jzAu + tfuv - b 1U + ^uW, 
d f v = — -uv. 

Furthermore, for A C d(c~ 1 D), we have uf^ D \A) = (a / c d )uff (cA) . 

The equation for (u,v) follows immediately by scaling the martingale problem (|2.1|2.2j) . The scaling of the 
exit measures follows by matching the finite variation parts of the semimartingale decomposition of Ut{4>) and of 
(a/c d )ubt(<p(-/c)) for <f> € C$(c D) given in the corresponding extended martingale problems. 



2.2 Decomposition results for the total occupation and exit measures 

Notation. For a solution (u, v) to (jl.ip we write uu it ] for the occupation measure f*u r dr (so that m s< t\(dx) = 
u(s,t,x)dx P a.s.) and write 

u \o oo) — hm ttrot], u(0,oo,x) = lim u(0,t, x), and u^f = lim uf D 

for the total occupation measure on D and its density, and for the total exit measure on dD. 

The lemmas below are stated under the restriction that D is bounded, which implies (see section [33)1 that the 
total occupation and exit measures are almost surely finite. The first lemma describes the monotonicity of the total 
occupation and exit measures with respect to the initial conditions (/x, /). We write = for equality in distribution. 

Lemma 4. Fix a bounded domain D and initial conditions fi = [±~ + /j + and f = f~ + / + . The law of the total 
occupation and exit measures of a solution (u, v) to on D with initial conditions (fx, f) can be decomposed as 

(«[0,=c),«£P) = (>,oo)+4,oo)'^+^ D ' + ) ( 2 -5) 

in either of the following two ways: 

(i) (u~,v~) is a solution with initial conditions (/i - ,/) and, conditional on cr{u~}, the process (u + ,v + ) is a 
solution with initial conditions = /i + and Vq = (where = hm^oo v^); 

(ii) (u~,v~) is a solution with initial conditions (p, f~) and, conditional on o~{u~}, the process (u + ,v + ) is a 
solution with initial conditions Uq — (3f + (l — cxp(—u~(0,oo)))dx and Vq = f cxp(— u~ (0, oo)). 

The next lemma shows monotonicity in the parameters j3, 7. 
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Lemma 5. Fix a bounded domain D, and parameter values f3 = (3~ + j3 + and 7~ = 7 + 7 + . The law of the total 
occupation and exit measures of a solution (u, v) to il.l\) on a bounded domain D with initial conditions (fx, f) and 
parameter values 7) can be decomposed as in \2. 5)) in either of the following two ways: 

(i) (it - , v~) is a solution with initial conditions /) and with parameter values f3~ , 7 and, conditional on o~{u~}, 
the process (u + ,v + ) is a solution with initial conditions Uq = /3 + (/ — v^dx and — and with parameter 
values (3, 7; 

(it) (u~,v~) is a solution to i f Li)) with initial conditions (p, f) and parameter values /3,7~ and, conditional on 
er{it~, v~} , the process (u + , v + ) is a solution with initial conditions Uq — 7 + U[g ^ and Vq — and with parameter 
values (3, 7. 

A special case of the Lemma [3] (i) will be particularly useful. The total exit measure can be built up in two 
stages: first run a process with p = (which is the well understood Dawson- Watanabe process) and then run a 
second solution started with the mass that the nutrient would have produced from the first process. 

The following simpler comparison will often be useful. It states, roughly, that if we convert some of the nutrient 
available into the equivalent amount of initial mass at time zero then we obtain a healthier process, in the sense of 
stochastic ordering. 

Lemma 6. Fix a bounded domain D and initial conditions fi and f,g. Let (u,v) and (ii,v) be solutions to 
11. 1]) on D with initial conditions (/i, / + g) and (p + f3gdx,f) respectively. Then for any bounded measurable 
F : A4(D) x M(dD) — > R that is non- decreasing in both variables 

E[F(u [0tOo)) u d J>)} < E[F(u l0tOo)) u d J>)}. 

The final lemma is a spatial Markov property, analogous to that known for Dawson- Watanabe processes. We 
choose two domains D~ C D + . We allow the boundaries dD~ and dD + to intersect, but we need this intersection 
not to be too complicated, and for this we ask that 

(dD~ n dD + ) n (dD~ \ dD + ) has surface measure zero in dD~ . (2-6) 

This is satisfied for example if both domains are boxes. We write [i\a for the restriction of a measure fx to a set A. 

Lemma 7. Fix bounded domains D~ C D + satisfying k2.6\) . Suppose (u,v) is a solution to 11. 1]) on D + started 
at (/x, /). Then the law of the total occupation and exit measures satisfies 

(w[0,oo),t4? + ) = («p,oo) + U [0,oo)> U '^ ) '~\dD+ + ) 

where (u~,v~) is a solution to hl.l)) on D~ with initial conditions (fiI(D~), fI(D~)) and, conditional on o~{u~}, 
the process (u + ,v + ) is a solution on D + with initial conditions 

u+ = u 9 J>->-\d++vI(D + \D-) and v+ = v~I(D-) + fI(D+ \ D~). 

Remark. As indicated in the introduction, these lemmas will be proved for discretized approximations to 
solutions, where the nutrient is used up in finitely many steps, and then via a passage to the limit. We spent 
considerable energy trying to find a proof using the continuum equation (jl.ip directly. Approaches using a change 
of measure (since the lemmas are easy to prove for a Dawson- Watanabe process) typically failed because the 
decomposition (|2.5p holds only for the total occupation and exit measures and fails for these measures at times 
t < 00. Another approach, using Laplace functionals, is sketched in remark (b) of section [3.11 On the other 
hand, we found the discretized approximations, in particular particle pictures, were the easiest testing ground for 
exploring possible extensions of these results to other models. 
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3 Change of measure arguments 



In section ETTl we review some results on Dawson- Watanabe processes and indicate when these imply estimates for 
solutions of (jl.ip . Section 13.21 gives the proof of Theorem [5J and section 13.31 contains the proof of the existence of 
the critical curve ^f(/3). 



3.1 Results from Dawson- Watanabe processes 

A Dawson- Watanabe process, with underlying Brownian spatial motion, killed at the exit of a domain D, and with 
a (deterministic) bounded, measurable mass annihilation/creation rate r\ : D — > R, is the solution on D, in the 
same sense as in section |2~T| to the equation 

d t u = Au - i]u + sfuW. (3.1) 

We shall abbreviate Dawson- Watanabe processes as DW processes. In the case r](x) = 7 € R is constant we will 
call the process a DW(£>, 7) process. We describe below some results for DW processes that we shall need, and we 
mostly indicate where to find the proofs, or how to adapt them, from the literature. The surveys by Perkins [TO] 
and Dawson [2] are our main sources. However, we give proofs for some of the results on exit measures, that are 
perhaps not in the literature, at the end of this subsection. 

A number of the estimates (for example moment bounds, estimates on the speed of the support, estimates on 
the extinction probability) also hold for solutions to (ll.ip . One way to see this is to set up a pathwise coupling for 
(u, v) a solution to (jl.ip on D so that ttj: < Ut < uf^ where «W is a DW(D, 7) process and is a DW(Z3, 7 — (3) 
process. We do not go through the steps to construct this coupling since, as we indicate below, the proofs of the 
desired estimates are typically based on stochastic calculus and are easily adapted to follow from the martingale 
problem for solutions to (jl.ip . 

1. Martingale problems. We quickly repeat here the martingale problem definition for solutions. A continuous, 
adapted Ai(D) valued process u, on some filtered space (f2, (J 7 *), P), is a solution to (|3.ip if for <f> 6 Cq(D) 

ut{4>) = u (cj)) + u s (A<p-r)<f>)ds + mt((j)), P a.s. (3.2) 
Jo 

where m t (<j)) is a continuous (JFt) local martingale with quadratic variation [m(</>)] t = J Q u s (<j) 2 )ds almost surely. 
Solutions, starting from /1 G A4(D), exist and are unique in law on C([0, 00), M{D)). Furthermore, for any solution 
there will exist an adapted process (uf D ) with non-decreasing continuous paths in A4(dD), called the exit measure 
process, so that for e C%(D) 

u t {4>) = u ((p) + f u s (A0-#)ds-uf D (0)+m t (0), Pa.s. (3.3) 
Jo 

where mt(4>) is a continuous local martingale with variation as before. The almost sure linearity mt((f> + ip) = 
mt((j)) +mt(ip), implies that if (f> n S Cq(D) converge bounded pointwise to ^ S C^(D) then m t ((f>n) —> m t {(j)) in 
probability (since [m(</> — 4> n )]t — * 0)- Now (|3.2p and (|3.3p imply that uf D ((f>) is measurable function of the path 
(u s : s <t). Choose is a countable family (4>k)k=i,... in C%(D) so that the map /1 — * {^(4>k))k=i,... is a continuous 
injection from Ai(dD) into R°° (with the topology of uniform convergence on compacts). Using the measurable 
inverse of this map, one one can construct a measurable map Rt : C*([0, t],A4(D)) — > Ai(dD) so that 

uf D = R t ({u s : s<t}) P a.s. (3.4) 

These maps show that the law of (u, u dD ) on is determined by the law of u on C([0, 00), M (D)). Also, for any 
DW process on D, the exit measure process can be constructed by using the maps Rt and then regularizing the 
path over the rationals. 
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Notation. We denote by Q® ,J] the law on of a DW process and its exit measure (u,u 9D ) on D, with mass 
annihilation/creation rate r\ as above, and starting at /i G .M(-D). We write Q?' 7 in the case that f]{x) = 7. 



The map /i — > Q^ ,V (B) is measurable for any Borel £? C Sl£>; furthermore the family of laws (Q®' v : /i G .M(-D)) 
form a strong Markov family in that for any solution u defined on (Q, (J 7 t ),J', P), and any finite (Tt) stopping time 
t, 

E[h(u T+ .)\T T ]=Q°?[h(U.)] Pa.s, 

for all bounded measurable h : C([0, 00), jM(-D)) — > -R. The measurability and strong Markov property follow from 
the Laplace functional (|3.6[) below by, for example, the arguments used in the proof of Theorem II. 5.1 of [TU] , 



For a DW process u, one can extend the local martingales to a local martingale measure, which we write as 
mt(<t>) = L J D (f> s ni(dx, ds), and is defined for <j> : [0, 00) x O x D — > R that are V x B(-D) measurable and satisfy 
L u s (4> 2 s )ds < 00 for all t > 0. The integral produces a continuous (.Ft) local martingale with quadratic variation 
[m((j))\t = /g u s ((j>1)ds. This follows as in [10] Proposition II. 5. 4. The same construction applies to give a local 
martingale measure for solutions to (|1.1|) . and this will be needed in our change of measure arguments. One simple 
use of these extensions is to allow time dependent test functions in the martingale problem. Since the exit measures 
uf D are non-decreasing, they induce a unique measure duf D (dx) on [0, 00) x dD characterized by the measure of 
the rectangles (s, t] x A being uf D (A) — u d s D (A). We will write the integral of ip : [0, t] x 3D — > R with respect to 
this measure as /„ du d s D (il>s). Then, if : [0, T] x D -> R is in the class C fc 1,2 ([0, T]xD), where and one partial 
derivative in time (denoted </>) and two partial derivatives in space exist and have continuous bounded extensions 
to [0, T] x D, the following decomposition holds: 



u t (<j)t) = u (M + J u s (A(/) s + <p s - ?70 s ) ds - J du d s D (cj) s ) + m t (4>) Pa.s. (3.5) 

This follows using easier versions of the arguments in [10 Proposition II. 5. 7. Moreover, the analogous extension 
works for solutions to (|1.1[) . 



2. Laplace functionals. Suppose first that D has a smooth boundary and that 77 : D — > R is smooth and 
bounded. The Laplace functional of the solution, the occupation measure and the exit measure is given, for smooth 
bounded h& : D -> R, : [0, t] x D ->■ R and : [0, t] x dD -> R by 



)A>7 



-u t (.hW)-j*u.{h^)-Siduf D (.hW) 



-M(0t) 



(3.6) 



where : s G [0,t], x £ D) is the unique solution, smooth on [0,t] x £>, to 



d 8 <t> 8 = A(/) s - ri4> s ~ + h 

M 



(2) 



on £>, 



/i^ s on dD for s G [0, t], and <fo = on D. 



(3.7) 



This follows by checking that s — > exp(— U s ((f>t- 8 ) — / ' U r (h s 2 "')dr — J Q * dU® D (h^)) is a martingale for s G [0,i]. 
If in addition 77 > and I? is bounded, the total occupation {/[o,oo) = hmj-foo J/[o,{] and total exit measures 
— lim t |oo 11®° arc finite almost surely and have Laplace functional 



-^[0,oo)(/n)-^ I '(fc2) 



-M(0) 



(3.8) 



for smooth bounded hi : £> — > [0, 00) and /12 : 91? — ► [0,oo), where (</>(x) : x G I?) is the unique solution, smooth 
on D, to 



A0 = — + rj(f) — hi on D, 



ho on 913. 



(3.9) 



This follows by examining the martingale t — > exp(—U t ((f>) — f*U s (hi) — U^ D (h 2 )). On bounded domains, and 
when 77 > 0, there exist finite solutions to (|3.9[) when hi,h 2 are negative and sufficiently small. In particular, the 
exponential moment Q^' v [exp(+9U^ ) (1) + 0U[o <oo )(l))] is finite for sufficiently small 9 > 0. 
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When the domain D is a box the associated p.d.e.'s have unique solutions that are smooth inside D and 
continuous on D, and the formulae for the Laplace functionals still hold, and can be established by an approximation 
argument. Indeed, one concrete way to do this is to suppose, without loss, that D contains the origin and to replace 
the test functions above by t (x/(l + e)) and <j>(x/(l + e), which are smooth in D. 



Remark (a). From the pde viewpoint, it is natural to use (|3.3p as a means of characterizing the exit measures 
(dx), that is purely as exit fluxes (as is suggested in jllj). Although we do not need such a characterization, 



t 



,dD 



since we have other constructions of the exit measure, we briefly sketch such an approach, and for convenience 
we consider deterministic initial condition /i and smooth bounded domain D. Construct the martingale measure 
m{dx,ds) from a DW(D, 7) process. Then for smooth h : dD — > R we may uniquely find <fi h £ C^(D) solving the 
Dirichlct problem Acf> h — on D and <p h = h on dD. Use the formula (|3.3p with the test function 4> h to define a 
continuous path t 



.dD 



(h) (although we have yet to establish that the values uf (h) arise from a measure u' 



,dD 



on dD). By adding the decompositions (13. 3p for cf) 9 and cf> h we find that uf D (g + h) 



.dD/ 



(g) + u% n (h) almost 



surely. Since ||<^||oo < Halloo, we can deduce a first moment bound sup s<t E{\u® D (h)\] < C(ju, t, 'y)||ft.||oo from the 
corresponding bounds for u. For <j) £ C%(D) with <j) = h on dD, the decomposition (|3.3|) holds by adding the 
decompositions for <fi h and <f)— (jr. Then the first moment bounds allow the decompositions to be extended to time 
dependent test functions <fi(t, x) as above. In particular we find, for h>0, that E[exp(-Xuf D (h))} < exp(-/z(<£ t )) 
where d t 4> = A0 — jcf) — on D, and <j> t = Xh on dD, and </> = 0. The inequality here arises since we do not yet 
know that uf D (h) > and so we have to localize via the stopping times r„ = inf{t : u t (l) > n}, and we obtain an 
upper bound by passing to the limit using Fatou's Lemma. However letting A — + 00 now implies that uf D (h) > 
almost surely. It is then not hard to establish the existence of a measure uf D so that u® D {h) = J dD h duf D . Finally 
the Markov property implies that u° 



dD 



,dD 



is non-negative and hence the paths of uf D are non-decreasing 



Remark (b). Mimicking the calculus that leads to the Laplace functional above, one can show that for solutions to 
(HUD that 

cxp \ -u t {4>) - U[ ,t](/ii) - uf D (h 2 ) + 13 J u s (v s (f>)ds\ 
is a martingale when <fr solves (I3.9[) with 77 = 7. One deduces, by letting t — > 00, that 



£ e -«[O.oo)('l2)-«f o D ('l3)+/3(0/4-C X p(-t l (O, ! X)))) 



-m(0) 



We do not know if this formula, as one varies h\,hi, characterizes the law of (tt[o )00 ),ti^). If it did characterize 



the law, a simple proof of all the decomposition lemmas from section [2721 would result, since it is straightforward 
to use calculus for both parts of the decompositions to yield the above formula for the sum. 



3. Some moments. Estimates on moments for the total mass process m*(1), starting from a deterministic initial 
condition, follow by choosing (j) = 1 in 



using uf (1) > 0, by the usual arguments (localizing by the stopping 
times t„ = inf{t : > n} and applying Gronwall's inequality). By rewriting (|3.3p , again with <p = 1, as an 

equation for uf D (l), one can then deduce moments for the exit measures. These arguments show, for instance, 
that for for p > 1 



sup0/ t (l)) p 

t<T 



U$ D {1)\ P 



<c( P ,r, 7 )(( / i(i) + ( M (i)) p ). 



(3.10) 



In particular, starting from a deterministic initial condition the processes mt(<fi) are true martingales. The same 
arguments apply to solutions to (jl.ip . and lead to 



S up(u t (i)y + \u° D (i)\p 

t<T 



< Q^~ 



sup(0i(l)) p + If/fW 

t<T 



(3.11) 



(We abuse logical order here, and below, by writing Q^'t' 1 before checking that these laws are uniquely defined. 



In such cases we mean that the arguments apply to all solutions to (II. ip .) 
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We will need the first and second moment formulae, for measurable g, h > 0, 

rt/\s 

Q D ^[U t (g)U s {h)\ = KGf n 9MG?' J h)+ rtG?«(G?«gG?«h))dr, (3.12) 

Jo 

where G^ > ' 1 g(x) = f D G® n ' x (y)g(y)dy and G^ n ' x (y) is the killed Green's function Gf n ' x (y) = e~ 7 *Gf ' x (y) for 
the domain D (that is Gf' x (y) is the transition density of a Brownian motion killed on its exit from D). These can 
be proved by using the time dependent test functions and $- h \ where <jy s = G^Slg for s G [0,t], in (|3.5p . For 
smooth g, h compactly supported in D, this test function has the required regularity, and the moment formulae 
can then be extended to general g, h by monotone class methods. In particular the same methods apply to the 
martingale problem for solutions to (jl.ip and give the bounds, for g > and k = 1,2, 

Q° n [(Ut(g)) k ] < Q%f [(U t (g)) k ] < Q^~ [(U t (g)) k ] , (3.13) 

illustrating the natural intuition that less killing leads to larger solutions. 

4. Extinction probabilities and speed of propagation. The probability of extinction for DW(i? d , 7) processes 
is given by Q R ,7 [Ut — 0] = exp(— A| 7 ' ) /i(l)) where Ai 7 ^ is determined for s G (0, t] by 

A (7) = -i(A (7) ) 2 - 7 A (7) , and Ai 7) -> 00 as s | 0. 

This follows by examining the martingale s — » exp(— Aj_ s f/ S (l)) for s G (0, t] (see the derivation of (|3 . 15(1 below). 
This same argument applies to solutions of (jl.ip to give the bounds 

Of a ~ P \Ut = 0] < Q°J n [U t = 0] and Q?/«[U t = 0] < Qf < 7 [t/ t - 0]. (3.14) 



DW processes, started from finite measures, are compactly supported at all times t > 0. The same will apply 
for solutions u to (jl.ip . by the absolute continuity of the law of u with respect to a DW process (see section l3~2|) . 
However we need a quantitative bound on the size of the support and we use two estimates for this behaviour, 
based on the results from Dawson, Iscoe and Perkins [3J. For R > let Dr = (—R,R) d . Choose < £r < 1 
smooth and satisfying {£r > 0} = IT^. Then for M > and smooth bounded 77 : R — » R, let S = (f>l M ' R ' r '' be the 
unique non-negative solution to 

= - - ^o = M^. 
< 00 as M | 00 an< i moreover <j)\ s R,ri ^ (x) — > as i? — > 00 (use a comparsion 



Then, for s > 0, 



1 



argument and [3] Lemma 3.6). Taking expectations of the martingale exp(— U s ((f>[^fg R )) at times s = and s = t, 
and then letting M — > 00, one finds that Q^ ,ri [Ut(D c R ) = 0] = cxp(— ii{4>\ R ^)). The same proof shows, when 
7 — /3/ > ry, the bound 



(3.15) 



Letting R — > 00 shows that solutions to (jl.ip started from compactly supported /i G A^(i? d ) must have compact 
support at any time t > 0. 



Similarly, let tp s 



^ e unique non-negative solution to 
ip = Aip — i-0 2 — + M£r, ip = 0. 



Then, for |af| < i?, %pi M,R ' 71 \x) f tpi Rv ^ (x) < 00 as M T 00 and moreover il)^ ,R, ^\x) — > as i? — > 00 (argue as in 
[3] Theorem 3.3 (ii)). Examining the process exp(— U s (ip^fg R )) as above one finds that Qf ,7? [J7[o,t](^D^) = 0] = 
cxp(— ^(^t )) an d the same proof shows, when 7 — /?/ > 77, the bound 



^[o,*]TO = ol > Qf [tr [0 , t] (S^) = 0' 



(3.16) 



11 



Letting R — > oo shows that solutions to started from compactly supported fj, remain compactly supported at 
all times. Comparisons can be used to estimate ip^ R ' r, \ For example suppose that Aw < (l/2)u> 2 — i]w on Dr and 



mi{w(x) : x £ dD R _ 
inside Da 



} | oo as e J, 0. Then il)\ R <won Dp for all t so that, when 77 > f3f 



C^[o,oo)(-Di 



> Q 



R ,n 



^[0,00) (-D j 



> exp(-/x(w)). 



7, for /1 supported 
(3.17) 



5. Occupation density. Define a putative density U(s,t,x) for C/[ s ,t] on path space £Id by setting, for < s < 
t, x £ D, 



U(s, t, x) = liminf J U r ([x, x + e) d D D) dr. (3.18) 
We also set ?7(0, t, x) = lim s | U{s, t, x). 



Sugitani [13] studied the occupation measure under the law Q R »° of a DW process on R d . His main result 
extends to our processes on more general domains, and with bounded annihilation/creation rates, as follows: Qa' v 
almost surely, the process (s,i,x) — * U{s,t,x) is continuous over < s < t, x 6 D, and x U(s,t,x) acts 
as a density for Uua. One way of representing the density is via a Green's function formula, as follows. Let 

G [o% = Jo Gf' x ds and G°'Jf = J* G^' x ds. Then, for a fixed t > S > and x G D, 



U(5,t,x) = U 5 (G°'*_ s] )+ I j^G°l*_ s] (z)(-ri(z)U s (dz)ds + M(dz,ds)) Q^" a.s. (3.19) 



where M(dz, ds) is the martingale measure associated to U . Moreover, for k € N, 

Q°' v [U(6,t,x) k ] < C(k, 6,n,\\ tjHoo.T) for all x € D,t € [5, T}. (3.20) 

When the initial condition is more one may take S = in ()3.19p . and indeed this is the case if \x has a bounded 
density. In particular, when /i(dx) = pdx we have, for k S N, 

Q%d V x [U(0,t,x) k ]<C(k,P,\\v\\oo,T) for all D,tE [0,T]. (3.21) 

It is not hard to extend Sugitani's proof of the representation formula to hold for domains D and with annihi- 
lation/creation terms, and also to bound the moments (inductively in k). Sugitani's proof of continuity is via 
Kolmogorov's continuity criterion and moments of increments. We omit the details since we will follow Sugitani's 
arguments to establish a similar increment estimates for approximations to solutions in section [4~2l and some of 
the key steps are illustrated there. 

The second moment formula (|3.12[) lead to the moment 

Q» n [U 2 (0,t,x)} = (MG&l")) 2 + /V(Gf' 7 ((G^]) 2 ))^- 
We use the bound Gf n ' x < e~ 7 *Gf ' x and the estimate 

q(t) := J (G^;*(y)) 2 dx < C*(d)i (4 " d)/2 for all x, y when d < 3. 
One can now deduce that 

ft 



Q 



D,7 



(U(0,t,x)) 2 dx 



D 



< C(d, 7) ^ 2 (l)q(t) + /i(l) / q(s)ds . (3.22) 



Again, one has that Q»f [(t/ 2 (0,<), 1)] < Q^-? [(U 2 (0,t), 1)]. 

6. A construction of the exit measure process and the spatial Markov property. A convenient way to 
construct the exit measure processes was developed in Perkins [TT] , which will be useful in our proof of a coupling 
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version of the spatial Markov property lemma below. This exploits a historical DW process, which is a process 
of random measures on the space of paths, where intuitively the paths trace the positions of the ancestors of the 
particles in a DW process. To explain the construction, we recall some notation from [TT] for the historical Brownian 
DW process (which is only required in this one subsection): on the path space C([0, oo), R d ) let C be the Borel 
(T-field and (Ct) the canonical filtration; A4(C) is the space of finite measures on C([0, oo), R d ) with the topology of 
weak convergence; on = C([0, oo), A4(C)) let Tt be the Borel sets and (H t ) the canonical coordinate variables. 
Write Q\ „ for the law on (f2#,7i) of a historical Brownian DW process, with branching rate one, with bounded 
measurable annihilation/creation rate r\ : R d — > R, and with initial condition /j € Mp(R d ). Letting AT be the Q V H 
null sets we define a filtration by Tit = D s>t (j{H r : r < s} V Af. The canonical process H = (H t ) is a historical 
DW process with creation/annihilation 77 on (f2# , H., H. t , Q% S)- In the notation from [TT| , it solves the martingale 

problem (MPf{% Ififj , where T)(s,y) = rj(y s ). 

For (Ct) predictable C : [0,oo) x C([0, 00), R d ) — > R, having left continuous paths of bounded variation, 
Theorem 2.23 of Perkins [TT] constructs the integral of C along the paths of a historical DW process as the limit 
(in probability, or almost surely along a subsequence) 

/ H s {dC s ) = lim V/(i" < t)H t n(C t n - C t » ), where tf = i2 n . (3.23) 

Note that the process J* H s (dC s ) is non-decreasing whenever C s is also non-decreasing. 

Fix a bounded domain D and define r = r(y) = mi{t > : yt ^ -D} for y £ C([0, 00), i? d ), and where 
inf{0} = 0. Define a measure u t (dx) on I? by 



u t (dx) = / 7(j/ t e dx, t < t, t > 0)H t (dy). 



IC([0,oo),R d ) 

If \x is supported inside D we may omit the indicator J(r > 0). For bounded measurable <f> : dD — > R define 

u a t D ($) = [ H a (dCt), where Cf(y) = <j>{y T )I( S > r > 0). 
Jo 

By definition u dD ((f)) > if <j> > and u dD (4> + VO = u dD (4>) + u dD (ip) almost surely. Using these properties, it is 
not hard to show that there is an (Tit measurable) random measure uf D on dD so that uf D ((j>) = J dD <j)(x)du d t D (dx) 
almost surely. 

Let Z(dy, dt) denote the martingale measure associated to H (constructed for example in section 2 of [TT]) and 
use the notation Z°(dy, ds) = Z(dy, ds) — r](y s )H s (dy)ds, that is an integral against Z°(dy, ds) means the difference 
of two integrals. This notation is borrowed from [11] and will save space since these two integrals frequently appear 
with the same integrand. We may apply the historical Ito formula ([TT] Theorem 2.14) to the process Y t (y) = yt/\ T , 
since the stopped path ytAr- For cj> e C^' 2 ([0, 00) x R d ) this yields 

4>t/\T(yt/\r)H t (dy) = fj,(cj) ) + J J ((p + A(p) sAT (y sAT )I(s < r)H s (dy)ds 

t 



4> S Ar(y S Ar)Z (dy,ds) fori>0, a.s. (3.24) 
Using Proposition 2.7 from [TTJ we have also, Q\ a.s. for all t > 0, 

Myo)I(r = 0)H t (dy) = fx(I(D c )<p ) + J J Mvo)I(r = 0)Z°(dy, ds). (3.25) 
Finally, Theorem 2.23 of [UJ shows, when C?(y) = <j) T (y T )I(t > t > 0), 

H t (Cf)= [ H s (dCf)+ [ [ (f> T (y T )I(s > t > 0)Z°(dy,ds) for t > a.s. (3.26) 
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0)H t (dy)~H t (Cf). 



u t {4> t ) = n(faI{D))- f H s (dCf)+ f [</> a (y a )I(8<T t T>0)Z(dy,ds) 

Jo Jo J 

+ J J \<j> + A(t>-r}(j)) s {y s )I{s <r)H s (dy)ds for t > a.s. (3.27) 

An approximation argument shows that (|3.27|) continues to hold for <f> £ C 1,2 ([0, t] x D). Taking time independent 
<j> shows that the extended martingale problem (|3.3|) holds for (u,u dD ). To check the regularity of the paths, note 
that the path t — > ut(4>) is continuous for <f> £ Cq(D) by (|3.2[) . To establish continuity when (f> £ C%(D), we use an 
estimate of the amount of mass near the boundary of D. For e > let (j>r> e (x) = (l — e~ 1 d(x, dD))\/0. If n(dD) = 
then U t (dD) = for all t>0,Q* * almost surely. This follows since t — > U t (dD) is continuous (see [TU] Theorem 

III. 5.1 and use absolute continuity between Qjf i7? and Q^' ) an d has zero first moment. Since 4>D,e I I(dD), for 
such \i we have 

sup t<T U t (4>D,e) | as e | 0, Qjj ,7? almost surely. (3.28) 

Note the simple coupling 

ut < u t for all * > (3.29) 

where u is a DW(i? d ) process with creation/annihilation rj, both with the same initial condition. This is natural 
since u does not have the killing on the boundary of D, and follows by setting u t {dx) = J I(y t £ dx)H t {dy). This 
coupling ensures that (|3.28|) is also true for sup t<T ut(0D,e)- This implies the continuity of t — > u t {<j>) for (f> £ Cjj(D), 
and hence of t — > u t £ M(D). The continuity of t — > uf D (h) for ft, € C^(dD) follows from the continuity of all 
other terms in (|3 . 3[) . This implies the paths t — > uf 13 are continuous, and they are non-decreasing by definition. 
So u is a DW process on D, with annihilation/creation rate r\, with exit measure u a£> , and with initial condition 
Hl(D). 

Dynkin [5] established a spatial Markov property for DW processes. We will need the following slight variant 
of the standard statement. The lemma is reasonably clear at the approximation level of discrete branching trees, 
but we choose to give a proof in the continuum setting using historical calculus. 

Lemma 8. Let D~ C D + be bounded domains. Fix fi £ A4(D + ). There exists a coupling of three processes: u a 
DW(D + , j) process with initial condition fi; u~ a DW(D~ , j) process with initial condition /iI(D~); and u + which, 

conditional on a{u~}, is a DW(D + ,j) process with initial condition fiI(D + \ D~) + and moreover 

these processes satisfy the splitting 

(m[o,oo),u^ I5+ ) = (>,oo) + u [0,o O y u ^ > almost surely. (3.30) 



Note that 

u t (^t) = J <t>tAr(ytAT)Ht(dy) - / '.',,(//(,)/(- = 
Combining this with (13. 2413. 2513. 261) shows that 



Proof. Define = inf{t : y t £ D^}. As above, we suppose the canonical process (H t ) is a historical DW process 
with constant annihilation rate 7 on (Qh , H, Ti-t, Q H „). Define 

ut(dx) = J I(y t £dx,t<T+) H t (dy), u^(dx) = J I(y t € dx, t < t~ , r~ > 0) H t (dy) 

u 9 t D \<P) = f H s {dCt + ), u d t D ~>-{<t>) = f H s (dCt--), 
Jo Jo 

where Cf ,+ (y) = 4>(y T +)I(s > t + ) for bounded measurable <f> : dD + — > R and Cf'~(y) — <p(y T -)I(s > t~ > 0) 
for bounded measurable <\> '■ dD~ — > R. The earlier arguments leading to (]3 . 27[) show that this defines processes 
(u,u dD+ ) and (u~,u dD >~) having laws Q^ +n and Q^jm-j as desired. 



14 



There exist random measures I on [0, oo) x D + and \ dD+ on [0, oo) x dD + satisfying, for bounded measurable 
[0, oo) x D+ -> R and rp : [0, oo) x dD+ -► R, 

poo p 

m = / U t -r-(yt)I(te(r-,T+})H t (dy)dt 



[°° H s (dCf) where Cf(y) = yj T+ _ T -(y T+ )I(s > r+ > r"). 
o 



We will show that conditional on a{u }, the law of (I, I dD+ ) is that of (Ut(dx)dt, dU 9D (dx)) under Q D 

I*' 

This allows us to define (u + , u 9D :+ ) with the required law so that 

/>oo />oo 

!(</>)= u+(<f> t )dt and I dD+ (i>) = dv° D+ >+y>.). 
Jo Jo 

The splitting (|3.30p for the total occupation measures holds since 

/>oc /> y* oo /> 

u [o,oo)( A ) = ^ y e A,t < r-,r- > 0)H t (dy)dt = J J I(y t eA,t< r-)H t (dy)dt 



Hl(D°)+u™ '-\ D +' 



while 



«p,oo)(4) = j[ / J (y f e ^ r+)H t (dy)dt, I([0,oo) x A) 
Note that for A C &D+ 



= jf J I(y t eA,te(T-,r + })H t (dy)dt. 



poo poo poo 

u d ° + {A) = H s (dC?>+), u™->-(A)= H s (dC?<-), I dD+ ([0,ac)xA) = H s (d6f) 

In In In 



where 



C?' + {y) = I(y T+ EA 7 s> t+), C?>-(y) = I(y T - eA, S > T -> 0), Cf(y) = I(y T+ eA,s>r+> T ~). 

Note also that Cf'~(y) = I{y T + € A,s > t + = t~ > 0) since A C 9Z? + . Then the splitting for the total exit 
measures holds since 

= + Cf>~ + I{y T+ £A,s>t+ = 0) 
and H t (T+ = 0) = for all t > 0, Q J H ^ a.s. 

To identify the conditional law of (I, \ dD+ ) we will use the Laplace functionals. Choose non-negative, 
continuous bounded h^'^ 1 : [0, oo) x Z? ± and h^'^ : [0, oo) x dD^ which are zero for t > T. Let (fr^ € 
C^ 2 ([0,T] x D ± ) n C([0,T] x Z>±) be the unique non-negative mild solution to 

-</>* = A0± - 7 0± - t^) 2 + /i (2)l± on D ± , 4>% = and = M 3 '^ on aD±. 
We extend these by letting <j>f = for t > T. It is enough to show 



e -«o(^) 



where Uq is defined as /i!(D + \ D 



,dD~ ,- 



\D+ 



This in turn will be implied by 



E 



= p -M(*o ) 



(3.31) 



(3.32) 



for all /i( 2 )^",/i( 3 )' 
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A comparison of the martingale problem for time dependent test functions (13. 5p with ()3.27j) shows that if 
(f> : [0,oo) x 0D+ -> R and C+' + (y) = <t> T+ (y T +)I{s > r+) then J* H a (dC*> + ) = f* du d s D+ ((f> s ) (and a similar 

representation for a du d t D ~ integral). Suppose first that (f>~ is in c£' 2 ([0, T] x D~). Then using (|3.27p on the 
domain D ~, one finds that defined by 

5 - = e -t*r tfn-J? «r(^ a) --)'i»-/o d«? D "'-(hi 3) --) 

i _ _ rT - / , (2). - rT t 8D~ ,— (3) , — \ 

is a martingale which satisfies = e~ J o "* * i at ~Jo au t t,' 1 * J anc j 







s 7 ^7 (y s )^"(s < r~ , r" > 0)Z(dy, ds) for i > a.s. (3.33) 



An approximation argument shows (|3.33[) continues to hold when (f> € C 1,2 ([0,T] x D ) n C([0, T] x D~). We 
now develop a similar representation for the Laplace functional of (I, I dD+ ). Set 



F t = ((t-T-) + -(t-T + ) + ),y tAT+ )= U I(ae(r-,T + ])da, / I(a < T + )dy(s) 
For G C b 1,2 ([0,oo) x i? d ) we have, using the historical Ito formula again, 
<f>(Y t )H t (dy) = ti<h)+ f I ' 4>(Y a )Z°(dy,ds) 



+ J J <P(Y s )I(se(T-,T + }) + A(f>(Y s )I( s <T + )H s (dy)d s . 
On {t > t + } we have Y t — (r + — r~,y T +) so that 

</>(Y t )I(t > r+)H t (dy) = [ <b T+ _ T -(y T+ )I(t>T+)H t (dy) 



= H t (Cf) + / MVt-W > r + = r-)H t (dy). (3.34) 

Using Theorem 2.23 of [11], 

H t (Cf) = f H s (dCf) + f I 4> T+ ^ T -{y T+ )I{s >t+> T-)Z°(dy,ds). 



On {t < t } we have Y t = (0, yt) so that 

(j>{Y t )I{t<T-)H t {dy) = I 4> (y tAT -)I(t <T~)H t (dy) 



f MVtAT-)H t (dy) + I MVr-)I(t > T-)H t {dy). (3.35) 

By the historical Ito formula again 

MvtAr-)H t (dy) = fj,(<p ) + J J <p (y sAT -)Z°(dy,ds)+ J A0o(»«)/(* < r-)H s {dy)ds. 
We may combine some terms from (|3.34p and (|3.35l) as follows, again using Theorem 2.23 of [TT] . 
, ,An r -)I(t > T-)H t (dy) - J MVr-W >t + = T-)H t (dy) 
4> I(D+)(y T -)I(t>r-)H t (dy) 



H t (Cf oHD+) '-) + j <p I(D+)(y )I(t >r-= Q)H t (dy) 

H s (dCf° I ( D+ ^-)+v((f>oI(D + \D-)) + I [ 4> I{D+)[y T -)I{t > r-)Z (dy,ds) 
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The last six displayed equations hold for all f > almost surely, and a little book-keeping combines them to yield 
<j> t _ T -(y t )I(t e(T-,T+])H t (dy) + J H s (dCt)-L(ci> I(D+\D-)) + J H^dCf* 1 ^-)) 

c/> s _ T -(y s )I{s e (T-,T+])Z(dy,ds) + I I (0 + A0- 7 0) s _ T -(y s )/( s <E (r' ,r+})H s (dy)ds. 



If we now suppose that (j> + A(j) — — (l/2)0 2 + h = for some h € Cb([0, oo) x R d ) then S/" defined by 
3+ = eM+f^oI(D + \D-)) + J^H s {dCt oI{D+) -)- J <l>t-T-(Vt)I(t e (r" ,T+])H t (dy)) 



exp( 

/o Jo 
is a non-negative local martingale and satisfies 



H s {dCt) - f* I h s _ T -(y s )I(s e (r- ,r+})H s (dy)ds) 



E+ = 1 - J J ~t<t>s-r- (y s )I(s € (r~ , r + ])Z(dy, ds) for t > a.s. (3.36) 

An approximation argument shows that (|3.36[) continues to hold with <j> — (f> + and h = h^' + . With this choice, 
and extending = for t > T, we let t — > oo to find that 



3f" exp 



(+ M +(^)-I(^ 2) ' + )-I 9D+ (^ (3) - + )). 
Combining (J3T33J) and (|3.36|) we have 

e +n+(0+)-i(h( 2 )-+)-i aD+ (^).+) e -/ -«r{k« w '-)*-/ -<"'-(*{ 3) '-) 

= e-»fc) - J™ J S t -S+0+ T .{y t )I{t e (r-,r + ])Z(dj/, dt) 

- jf y S+H t -^ t -(y t )J(f < t-,t- > 0)Z(dy,dt) (3.37) 

(noting that the cross variation of the two stochastic integrals is zero). Note that Sj < exp(+||<?!) + || 00 u ( |(l)). The 
quadratic variation of the first stochastic integral on the right hand side of (13.37[) is therefore bounded 

||0+|| 2 oo exp(+2||0+|| oo ( M (l)+^--(l))n [o , oo )(l) 

The finiteness of small positive exponential moments for the total exit measure it^f implies that the stochastic 

integral has mean zero if the norm ||<^ + ||oo is small, which in turn is satisfied if the norms ||/i*- 2 - ) ' + |Joo and ||/i^ 3 - > ' + ||oo 
are small enough. Similarly the other stochastic integral is mean zero. Taking expectations in (|3.37|) establishes 
the conditional Laplace functional (|3.31[) when ||/iA 2 )' + ||oo and ||/i^ 3 - ) ' + || o are sufficiently small. But this is sufficient 
to characterize the conditional law. □ 



3.2 Existence and uniqueness via change of measure 

Proof of Theorem [2l Many of the change of measure arguments we will use can be found in section IV. 1 of 
Perkins [TU]. Fix / 6 B(D, [0, 1]), and on the path space {QdMMuQu' ) use tnc path space density from (|3.18[) to 
set Vt(x) = f(x) exp(— U (0, t, x)), where we use the convention exp(— oo) = 0. Note that V is V x B{D) measurable 
and has non-increasing paths as required. Define the stochastic integral 

Mf 7,/ = / [ (0V s (x)-j)M(dx,ds) (3.38) 

J JD 
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where M(ds,dx) is the martingale measure on [0, oo) x D constructed from U under Q^ ,Q (as explained in sec- 
tion 3.1.1). A simple version of the arguments in [10] Theorem IV.1.6.(b) implies that the stochastic exponen- 
tial Et{M^- n ^) is a true martingale, so that there exists a measure Q^f n on the path space (£ln,U) satisfy- 
ing dQ^j' 7 /dQ^' a — £t(MP' 7 'f) on Ut- It is straightforward to check via Girsanov's theorem that the triple 

(Ut,Ut D ,Vt) solve, under Q^f n , the extended martingale problem (12. 3|) and (|2.4|) for solutions to (jl.ip . The 
regularity of U(s,t,x) will imply that (|2.2[) holds, namely for all bounded measurable (j> 

(V t ,<j>) = (/, 4>) - ! U s (V s ^)ds, Cffi a.s. (3.39) 
Jo 

This would be trivial, by the fundamental theorem of calculus, were U s {dx) to have continuous density in (s,x). 
We show at the end of this section that it follows in our more singular setting. This completes the proof of part (i) 
of Theorem [2] 

To show uniqueness of solutions, we take a solution (u,v) to (jl.ip on D, on a filtered space (fi, (Tt),T,P), 
with an arbitrary Tq measurable starting condition. Let m{dx,dt) be the local martingale measure extending the 
local martingales m t {4>)- We now reverse the change of measure argument. The local martingale defined by 



if ' lJ = S t (J J (7 - (3v s (x))m(dx, ds) 



is in fact a true martingale (argue as in 1 1 Oj Theorem IV. 1.6. (a), noting that the deterministic initial conditions 
there are not needed for this argument). Defining Q by dQ/dP = mf' 7 '^ on Tt, one finds that (u t ) solves the 
martingale problem for a DW(D, 0) process under Q, and so Q is determined on a{u s : s < t} by the law of Uq. By 
Sugitani's results listed in section f3.ll the occupation measures it[ s ,t] have densities (under P or Q), and moreover 
that there is a density (t, x) — » u(s, t, x) that is continuous on {(s, t, x) : < s < t, x S D}. Knowing this one shows 
that v is given by (see the argument at the end of this section) 

v t (x) = v {x)e- u( -°- s < x) for all t > 0, for almost all x, P a.s. (3.40) 

where u(0, s, x) = lim^o u(r, s, x). Now we may replace the change of measure by 

dQ/dP = £t(J J (t ~ (3v (x)e- u{0 ' s > x) ^ m(dx, ds)^j on T x 



since the two exponential martingales are equal upto a null set. Inverting one finds 

dP = E t (^J J (pv Q {x)e- u(0 ^' x) -7) rh{dx,ds)^jdQ on T t . 

where rh(dx, ds) is the local martingale measure determined by the local martingales found in the martingale 
problem of u under Q. The stochastic exponential is a measurable function of (vq,{u s : s < t}). Applying the 
Markov property of u under Q at time t — 0, we find that P is determined on a{u s : s < t} by the law of (uq, i>o)- 
Indeed, for each t > there is a measurable kernel pt(fJ-, /, B) for / e B(D, [0, 1]), fi € M(D) and Borel B C M{D) 
so that 



P[iH €dB}= / pt(p, /, B)P[u a e d/i, v € c(f]. (3.41) 
Jb Jm(d) 

Under Q, and therefore under P, there exists an exit measure process u dD satisfying uf D = Rt({u s : s < t}) almost 
surely. Another Girsanov calculation starting with the extended martingale problem under Q, shows that the exit 
measures satisfy the extended martingale problem for (jl.ip under P. This completes the proof of parts (ii),(iii) 
and (iv) of Theorem For part (v), note that the measurability of the maps 

/) -» Q°f n ( B ) follows from the 

analogous measurability of Q^'° and the fact that, for t > 0, the derivative dP/dQ®'° on Tt is a measurable map 
of (/, {u s : s < t}). The strong Markov property for the laws Q can be deduced either from the strong Markov 
property for the DW laws Q®' 1 , or more easily from the uniqueness of the one dimensional distributions for the 
martingale problem (use (|3.41|) and follow the argument in [TO] Theorem II. 5. 6). □ 
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Remark. It is straightforward to show that uniform integrability of the martingale t — > E^M^' 1 ' 1 ) under Qj} ,0 

for some pL 7^ (or that of t — * £t{M^'°' x ) under Q^"*' 7 ) is equivalent to certain death for the parameters (J3, 7). 
Sec the remark in section 177X1 



Proof of (|3.39p . Suppose [u s : < s < t) is a continuous path in M{R d ). Suppose also that (n(s,x) : s < 
t, x G R d ) is bounded and continuous and that J v r dr = n(s,x)dx as measures for s < t. Let T$ denote the heat 
semigroup with generator A, and write Tg for the dual semigroup acting on measures. For continuous compactly 
supported ip, we shall argue that 

(1 - e~ n ^)i;(x)dx = lim J (1 - e - Tsn ^' x) )ip(x)dx 



lim 

.510 



e- TsnM <ip(x)Tg*v s (dx)ds 
e- n< - s ' x ^(x)u s {dx)ds. 



Since the measures TgV a have densities, the paths s — > Tsn(s,x) are absolutely continuous for almost all x. So the 
second equality above follows by expanding s — > exp(— Tgnl-s, x)) as an integral of its derivative over [0, t] and then 
interchanging s and x integrals. The third equality holds since Tsn(s,x) — > n(s,x) uniformly over s < t and x in 
the support of tp, and since Tgi/ S — * v s weakly. 

The assumptions above hold, Q®'j'~ l almost surely, if we take v s (dx) — ip(x)U e + s (dx), for some e > and 
continuous ip, with compact support in D, and take the corresponding density n(s, x) = ip(x)U(e, e+s, x). Choosing 
continuous ip compactly supported in D and ip = 1 on the support of ip we obtain, almost surely, 



-U(e,s,x) 



ip(x)U s (dx)ds. 



To pass to the limit e I in this equation note that for t > e 



-U(e,t,x) 



-U(0,t,x) 



dx < U(0, e, x)dx 



[ U s {l)ds^0 ase|0, 
^0 



< IIV-II 



and similarly 

rt f 

- u ^ s ' x) ^{x)U s {dx)ds- 
U s {U(0,e))ds 
U(0,t, x)U(0,e,x)dx 
which converges in C 2 to zero, as e J, 0, by (|3.22p . This leads to 

e- u ^ x) iP(x)dx = f ^{x)dx 



e- u (°' s > x) ip(x)U s (dx)ds 



< 



- u ^°' s ' x ^(x)U s (dx)ds 



for continuous compactly supported ip. We may extend to general bounded measurable ip : D 
class arguments, and choosing ip = f<j> for tj> £ C®(D) shows that (|3 . 39[) holds. 



R by monotone 
□ 



Proof of (I3.40|) . First, we may use a density argument to ensure (vt, <fi) = (vo, 4>) — J Q u s {v s (p)ds for alH > 
and all continuous <f> with compact support in P a.s. Now we can argue pathwise. Fix a sample point w so 
that the above happens and so that also (r, s, x) — > u(r,s,x) is continuous over x S D,0 < r < s < t. Choose a 
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partition 0<s = to<ti<...<tN=t with m.ax.i(U — U-i) = S. Choosing 4>(x) = exp(— u(t n -i,t, x))ij)(x), for 
some continuous ip > with compact support A C D, we find 



A little rearrangement leads to 



u(t n -l,t) 



-u(t n ,t) 



1> 



vt^e-^-^iiS-El + El 



where E l n > are given by 



tn-1 



^ e -«(i„_ 1 ,t) (Us _ Vt j^ ds < ii^u^ ( u ( tn _ u tn )I{A), Vt^ - v tn ) 



and, using < e z — 1 — z < 2z 2 for z 6 [0, 1], 

E 2 n = (vt^e-^-^^^-^ - 1 - u{t n ^,t n )) < 2||V||oc (/(A), « 2 (<„-i, <„) A 1) . 



Now sum over the partition to find 



^) = Ue-^H)-Y. E n+H E l- 



(Vt, 



The error bounds above, and the fact that max{u(s', t' , x) : s < s' < t' < t,\s' — t'\ < S, x E A} — > as S — > 0, 
imply that J2 n E^ — ► as (5 — * 0. Thus = w s (x) exp(— w(s, i, a;)) for almost all x. The non-increasing paths 

t — > ^(x) ensures that i>t(a;) = exp(— u(s, t, x)) for all < s < t, for almost all x. Since Hui s ^qv s (x) = vo(x) 
for almost all x (use the non-increasing paths of v and (|2.2[) ) we reach the conclusion (|3.40[) . □ 

We conclude this section with the proof of local extinction of solutions. 
Lemma 9. Suppose that 7 > 0. Then, for compact A C R d , 

Qf/' 7 [tf[t,oo)(^) = 0] -» 1 as t > 00. 

Proof. Although not directly used, we first give a weaker fixed t extinction result whose proof is similar to the 
support results in section l3~3l Fix R > and let < ipR < 1 be smooth and satisfy {ipu > 0} = Br — {x : \x\ < R}. 
Let 4>[ M \x) be the solution, for x G R d , t > 0, of 



Then stochastic calculus shows that, under Q„J , 



d t <j> = - 7^ - ^- + (34>I{t < 1), <j>o = M^ R . 



so that, 



f e -Us(4-l)\ > -^ e - a "(^-i){7 s (y s 0W)J(s < t - l)ds + martingale increments 



Q«J^{U t (B R ) = 0} = lim Q^J 



M- 



-MUMr) 



U s {V s <j> t -s)I{s <t-l)dt 



(3.42) 



where (f> = liniM^oo 4>^ ■ As in the extinction results in section 3.1.4, a comparison argument shows that <f> t is 
bounded for all t > 0. Moreover, then ||</>t||oo — > as t — ► 00, so that the term exp(— fi{4>t)) —* 1- Similarly, for 
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fixed (s,uj) we have U s (V s 4>t-s) — ► as t — > oo. To justify the conclusion that j^ n \Ut{B R ) = 0] — > as i — > oo 
we will dominate the second term on the right hand side of (I3.42[) by 



<P*(x)f(x)(l-e- u ^'°°'^)dx 



where <f>*(x) = sup t>1 <pt(x), and where we are argue as in (|3.39p for the first equality. But, for 7 > and where T t 
is convolution with the heat kernel pt{x), 

4>*(x) < supe~ 7 *T t 0i(a;) 



< sup 



y|<i 



<pi(x + y)e lt Pt(y)dy\+ / </>i(a; + y)sup{e 7 *p t (s/)}dj/ 



< sup{</»i(z) : |z - x\ < 1} + C(j) / + y)e- lU2 ^dy 



The fact that (</>*, 1) < 00 can now be deduced from the boundedness and exponential decay of </>i, which follow 
from comparison arguments as in [3] Lemma 3.1. 

The stronger conclusion of the lemma follows by a similar argument using a slightly more complicated test 
function. Suppose that ip R satisfies in addition ip R = 1 on Br_i. Fix to,tx > and redefine <fr as the solution to 



d t (j> = Acj) - 70 - + Mi/> R I(t < h) + (3^ R+2 I{t <h + 1), cb = 0. 

Some calculus shows that, under Q R , 

d(exp(-C7 s (0 to+il _ s ) - M[fy DitD v a ]W>ii))) 

= exp(-C7 s (0 to+il _ s ) - MU [tQttoVs R )) (-0U s (V s ct>t o+tl - s ) + 0U s (ct>t o+tl - s ip R+ 2)I(s >t - l))dt 
up to martingale increments. Taking expectations we have 

t7 a (^ to+tl _ a y a ) - U s (<f> t0+tl - s ^R+2)I(s > t - l)ds . (3.43) 



> e -M(0t o +ti) _ pQ Rd >P''~l 



There is a unique non- negative solution (see the arguments of [7]) to 

- - 2 - 

= 70 + — — P4>ip R +2 on and 0(x) | 00 as |a;| J. R. 

A comparison argument shows that 4>t(x) < 0(aO for t < t\ + 1 and |a;| > R. Also 4>t 1 +t < A t for i > where 

A 2 

A s = — 7A S ^ + (3X s I(s < 1) and A s | 00 as s J, 0. 

Moroever, as in the fixed i argument, 

</>t 1+ i+ s < e-^T s <j) tl+1 < e~ JS T s ((j) A Ai) := S . 



Using these bounds in (|3.43p we obtain for to > 1 

e -MU ltoito+H] (ipR) 



e 



cr.(V;^ _,_i)d»+ / U s (V s <f>I(B R+1 ))ds 



to-l 



Now we let Af, ii — * 00 to obtain the same upper bound for Q^/^' 7 [^[t Q ,oo)(-Sfj) = As to ^ 00 we have At 
The conclusion of the lemma follows as in the fixed t result using the domination J °° U s (V s (p )ds < {<f> ,1) where 
4> = ma,x{<fiI(B R+1 ), sup s S }. The integrability (tfi , 1) < 00 follows as in the fixed t result using the exponential 
decay of the function <f>. □ 
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3.3 Existence of the critical curve \I/(/3) 



For a DW(D, 7) process the death time time r = inf{£ : U t — 0} is finite, and then f/ro,«0 = ^[o,r] an d = D 
are finite measures, Q® n almost surely. Using the change of measure from section 1^21 the same holds for the law 
Q°j n when D is a bounded domain. Indeed the change of measure martingale M 13 ' ^ from (|3.38[) satisfies 

[MP' '*]. = P 2 f I f{x)e- 2U ^ s ^U s {dx)ds 
' Jo JD 

by arguing as in the proof of (|3 .39[) . This is bounded independently of t, for a bounded domain D, so that £ t {M^' '^) 
is a uniformly integrable martingale, ensuring the above almost sure properties carry over to the law Q^'j 3 ' 7 . 

A similar argument leads to the following lemma, which shows that either certain death or possible life must 
occur, for each pair of parameter values (3, 7. 

Lemma 10. If ' {^ n [Ut 7^ for all t] is strictly positive for some /1 G A4(R d ) then it is strictly positive for all 
non-zero n S M.{R ). 



Proof. Step 1: We claim that the laws of solutions y' 3,7 and Q^fg^' 1 are mutually absolutely continuous on 
the entire sigma field Uoo = a{Ut : t > 0}, provided that f and g differ only on a compact set. To see this, consider 
the Radon-Nikodym derivative between the two laws, which is the stochastic exponential £(M) arising from the 
martingale Alt = P f Jnd(f — g)(x)e~ u i°- 3 ^ x " , M(dx,ds). Then, if < /, g < 1 are supported inside the compact 
set A, 



[M] t = f? [ [ (f-g) 2 (x)e~ 2U ^^U s (dx)ds 

Jo JR d 



< P 2 e- 2U i°^U s (dx)d 



s 



J A 



2 



(l - e- 2U l°^ x ^ dx < (3.44) 

The exponential martingale is therefore uniformly integrable which establishes the claim. 

Step 2: A result of Evans and Perkins [6 j shows that, for any fixed t > 0, the laws Q®'°[Ut € •] and 
Qn'°[Ut G ■] °f two DW(R d ,0) processes started at non-zero /1, fi are mutually absolutely continuous on M.{R d ). 
The change of measure ideas used to show existence and uniqueness imply that the same absolute continuity holds 
for Qff'-'pt G •] and Qff n [U t G ■]. 

Now we combine the two steps. Suppose [i, fj, are non-zero and compactly supported, and that started at 
(/i, 1) the solutions may survive. Then 

< Q£i'"' 7 [I7 t ^0faraM 

= / Qf g ^[Ut + for all t] Q«? n [U t0 G du, fe~ u ^ G dg] 

JB(R d ,[0,l])xM(R d ) 

by the Markov property at time t$ > 0. By the compact support property of solutions, / exp(— E/(0, to)) is identically 
one outside a compact set, almost surely. So, using step one, 

< / Qf/'^[U t ^OhTallt]Qf^[U to edu,fe- u ^edg] 

JB{R d ,[0,l])xM(R d ) 

Qfi'^Pt ¥= for all t] Qff n [U t0 G dv]. 

M(R d ) 
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Now using step two above we obtain < J M Q vl ' 7 [£/t 7^ for all t] Qn{ ^[Ut a S dv] and undoing the steps we 
find that the process with initial conditions (ft, 1) may survive. 

To remove the restrictions that //, ft be compactly supported, recall that the solutions are compactly supported 
at any time t > almost surely. Hence, if solutions may live from some general /i ^ 0, they may, by applying the 
Markov property at time t > 0, live from some, and hence all, non-zero compactly supported initial conditions. 
Any solution started at non-zero ft will have positive probability of being alive and compactly supported at small 
times t > (use the continuity of the total mass process) and hence, by the Markov property again, may live for 
all time. □ 

The next two lemmas give a characterization of certain death, which allows us to reduce the question of 
possible life/certain death to the study of the total exit measures solutions on bounded domains. The existence of 
a non-decreasing critical curve follows from this, since these exit measures are stochastically monotone in /3, 7. 

Lemma 11. Suppose /i € A4(R d ) is compactly supported. Then y' 9 ' 7 almost surely 

{Ut = for large t} — {J7[o,oo) is compactly supported}. 

Proof. The inclusion {Ut — for large t} C {J7[o,oo) is compactly supported} follows from the compact support 
property (|3.16[) . Let Dl = {—L, L) d . We claim for any L, P a.s., there exists a t > so that either Ut(dx) = or 
Ut{D c L ) > must occur. This implies the opposite inclusion. 

Applying the first and second moment bounds (|3.12I3.13|) . and the simple estimate P[Z > 0] > (E[Z)) 2 / E[Z 2 ) 
for a non-negative variable Z, we may find t\ > 0, depending only L,(3, 7, so that Q^'f n \Ui(D c L ) > 0] > ei for 
all / e B(D, [0, 1]) and > 1. Using the extinction probability (|3.14p we may find £2 > 0, depending only on 
/3,7, so that Q^'^' 7 [[/i = 0] > £2 for all / 6 B(D, [0, 1]) and < 1. Applying the Markov property at time t, we 
obtain 

Qf/ a [U t+ i(dx) = or U t+1 (Dl) > | U t \ > e 1 A e 2 . 
Iterating this estimate over the integers t = 1, 2, . . . and applying Borel-Cantelli completes the argument. □. 
Lemma 12. Suppose fi 6 Ai is compactly supported in D. Then 

0%"" [f/[o,oo)(C c ) = 0] = Q D J^[U d J> = 0]. 

Proof. We may construct a coupling of solutions (u~,v~) and (u + ,v + ) as follows: (u~,v~) is a solution to 
(jl.ip on D with initial conditions (/i, 1); u 9D '~ is the associated exit measure process; r = inf{£ : u d t D '~ > 0}; 
and conditional on o~{u~~ ,v~} the process (u + ,v + ) has the law of a solution on R d started at the random initial 
condition Uq = u~ and Vq = v~. One way to do this is to use the measurability of /) — * Q^f' 1 and construct 
a skew-product measure on the product space £Id X ^r^- We may also suppose that i> t _ = exp(— u~(0, t)) and 
w t + = cxp(— u~(0, t) — u + (0,t)) for all t > 0. Now define (u t ,v t ) — (u^,Vf) for t < r and (u t ,v t ) = (u^_ T ,v^_ T ) 
for t > t. Then w f is a measurable function of (u s : s < t) and this can be used to check that, with respect to 
the filtration a{u s : s < t}, the process v is predictable and non-increasing. For <p G C 2 (R d ), we may apply the 
extended martingale problem for (u~,v~) to the restriction of 4> on D. Combining with the martingale problem 
for (u + , v + ) we find that (u, v) is a solution on R d with initial conditions 1). Moreover on the set {r = 00} we 
have that W[o,oo)(-D c ) = u^ Q oc - ) (£ )C ) = 0. Thus 

Q*?< 7 [^[0,0.) (^ c ) = 0] > P[r = co] = Q D J ' 7 [D« = 0] . 

For the converse inclusion we construct v~) a solution to on i? d with initial conditions (/x, 1); r = inf{i : 
u^(D c ) > 0}; and conditional on er{u~,v~} a process (u + ,w + ) with the law of a solution on D started at the 
random initial condition Uq = u~ and Vq — v~. A result of Perkins (see [TU] Theorem III. 5.1 and use absolute 
continuity) implies that the paths of t — ► u^(D c ) are almost surely continuous. So we may define (u t ,Vt) as above 
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and the process u is almost surely continuous. A test function <j> € Cq(D) may be extended to <ft € C£(R d ). 
Applying the martingale problem for (u~,v~) to <p and the martingale problem for (u + 1 v + ) to <f>, we find that 
(u, u) is a solution on D with initial conditions ([i, 1). Moreover, choosing uf D = UnjX on t > r and itf D = for 
t < t one finds the extended martingale problem is solved. Then on the set {r = oo} we have that = and 
thus 

Q°f" = 0]>P[t = ™} = Qj> [C/ [0 ,oo)(^ c ) = o] 
completing the proof. □ 

Using the above results and the comparison results stated in section [2~2l we will now deduce the existence of 
an non-decreasing critical curve '5 as in Theorem [1] 

Corollary 13. There exists an non- decreasing function *5f((3) : [0, oo) — > [0, oo] so that for < 7 < ^((3) possible 
life occurs and for 7 > ^(/3) certain death occurs. 

Proof. Lemmas [TUl ED and QJ show, for /i ^ 0, that Q^f'^pt = for large t > 0] = 1 if and only if Q^f^pt = 
for large t > 0] = 1 if and only if Qf g ^' 7 [t/[o,oo) is compactly supported] = 1 if and only if 

supQ^[l^=0] = l. 

n 

But Lemma [5] shows that these probabilities are non-increasing in f3 and non-decreasing in 7. The result follows 
by setting * (/3) = sup{ 7 > : sup„ Q^ 7 [[/^" = 0] < 1}. □ 

Remark. The change of measure (|3.38[) can be used to show that (/3, 7 ) — > Q^f' 7 [U^P = 0] is continuous when 

Z? is bounded, and hence that Q^^ n \U t — for large f] is lower semicontinuous in /3, 7. This does not seem to 
have any immediate implications for the continuity of ^ . 

4 Approximations 

The proof of the decomposition results in section 12.21 is rather clear for particle systems that approximate our 
reaction diffusion system. Our first proofs used a full particle approximation with population and nutrient particles 
living on discrete lattices. We later realized that the key to establishing the comparisons was to discretise the effect 
of the nutrient, and we here present an approximation where the reaction with the nutrient occurs in a finite number 
of discrete 'packages' but the population still evolves as a continuum SPDE. This intermediate approximation makes 
passage to the limit easier to establish. This passage to the limit is broadly similar to many in the literature, and 
the points of most interest are perhaps (i) that the nutrient interaction is singular as a function of the population 
measure u, and requires tightness of the occupation densities; (ii) that convergence of the exit measures follows 
quite simply from convergence of the population measures, via the extended martingale problem (|2.4[) . 



4.1 Construction of the approximation 

The approximations will depend on a parameter N, which we will suppress in the notation in this section. Fix a 
domain D C R d and partition it as a disjoint union of sets D = UDj, where each set Dj has diameter at most 
A^ 1 . Choose a finite number of functions (ipk '■ 1 < k < K), each of the form N~ 1 I(x € Dj) for some Dj in the 
partition. Each ipk will represent a small package of nutrient which may be triggered to produce new population. 
The function / = ipk < 1 will be our approximation to the initial nutrient level. 
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Choose a probability space equipped with the following independent families: an i.i.d. family of rate one 
exponential variables &k for k > 1; independent DW(D,7) processes (uk,t '■ t > 0), for k > 0, with initial conditions 

u o,o — an d Wfe.o = P~>Pk(x)dx for fc > 1. 

Let (itf ^ : t > 0) be the associated exit measures processes. Given realizations of these variables, the approximation 
can be constructed pathwise. We will list as St C {1, . . . , K} the labels of those nutrient packages that have been 
triggered by time t, and denote by Tk € [0, oo] the time at which the fcth nutrient package is triggered. Thus So = 
and St — {k : Tk < t}. The approximation will satisfy 



ut = u 0it + u k ,t- Tk , 
kes t 



i)D 



keSt k£St 



(4.1) 



Moreover we want the approximation process to be triggered in such a way that 

T k = inf jt : U[ 0t t\(ipk) > efcj for k > 1 

where "0^ = ipk/(ipk,^-)- This uniquely specifies a process, which can be constructed pathwise by defining the 
triggering times and the processes u t ,v t inductively over the intervals between triggers (and where, on a null 
set, packages may be simultaneously triggered). 

We now derive a martingale problem for the approximation. On the intervals [Tfc,Tfc + i) the approximation 
process evolves as a finite sum of DW processes. So it satisfies the following martingale problem: for cf> E C 2 (D) 



= n{4>) 



f u s (A<f>~^)ds~uf D (((>)+ V^ sMs (0) 
Jo s <t 



"it (</>), 



(4.2) 



where mt((f>) is a continuous martingale with quadratic variation J* u s (<f> 2 )ds and where D s u s ((f>) is the jump 
u s (4>) — u s -((j)) at time s. We can compensate the jump term as 



s<t ^° fe05 s 



(4.3) 



where E^\<p) is a martingale with jumps bounded by j3N 

by 



-d-l\ 



and with predictable brackets process given 



d(£«(0)) t = (3 2 Y,(<t>,A) 2 uSk)dt 

k&St 

k<£St 



UWlcN-^u^dt. 



(4.4) 



The term N arises from bounding (1, ipk) < N (since < ip k < N 1 and the diameter of the support of 
ipk is at most iV _1 ). Examining the variation of the function <f> over the support of ipk, we can approximate 

(4>,ipk)u s (i(>k) - u s (ip k <P) < N^WV^WooUsiipk) 

and rewrite the compensator in (14. 3[) as (i J Q u s (v s (p)ds up to an error of size O^N^ 1 ). Combining these estimates 
we have 

(4.5) 



ut(<l>) = n(4>) + / u s (A<j 3 + l3v s <j 3 ~ 1 (j ) )d S ~u d t D {4>)+mt{(j})+E[ 1 \ ( t ) )+E[ 2 \4>) 
Jo 

(2) 

where the second error term E\ (</>) is controlled by 



£ t ( %) - E[ 2 \ct>) 



ks?S r 



< /? ((jV — 1 1| V0|| oo) A ||0||oo) / u r (l)dr 



(4.6) 
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Taking cj> = 1 in (14. 5|) we find that 

u t (l)</i(l) + (/3- 7 ) / u s (l)ds + m t (l) + £ t (1) (l). 



From this, via the standard argument (localizing, applying Burkholder's inequality to control the martingales and 
Gronwall's inequality), one can derive the moment bound, for any p > adnd T < oo, 



E 



su P Ml)P> 

t<T 



<C(p,/3, 7 ,T)(l + ( M (l)) p ). 



(4.7) 



We now consider the martingale decomposition of t — ► Vt (x) . For x G D let be the partition element that 
contains x, let ipx{y) — N^ 1 I(y G D x ) and ip x = ip x /(ip x ,^}- Then a nutrient package satisfies ipk(%) > only if 
ipk = ipx- The process t — > i>t(x) is a pure jump process with compensator given as in 



^0*0 = f( x )- i)k(x)u s {j)k)ds + m%(x) 



= f( x )~ / ^(x^^Jds + m£(x) 



for a martingale mj" (x) with jumps bounded by N 1 and with predictable brackets process 

(x)u t (ipk)dt < N v t (x)u t (ip x )dt. 

We can solve the equation (14. 8|) for v t (x) as 

Vt (x) = e -«[o.«i(tf.) C/(a;) + / e""i»-i (i, (imj(a:)) . 



(4.8) 



(4.9) 



(4-10) 



4.2 Tightness of the approximations 

To pass to the limit in the martingale problems for the approximation we require tightness estimates. We tried 
deriving these from known tightness estimates for each of the DW processes used as building blocks, however the 
random start times caused some difficulties in combining the estimates, so we proceed by repeating Sugitani's 
estimates for our discretized approximation. 

We now include the dependence on N and write the approximation constructed in section |4~T1 as (u^ N \ i)( N >). 

Proposition 14. Let D be a bounded domain. Suppose the approximation (u^ N \v^ N ^) has initial conditions 
u^ N) = and v { N) = / W < 1 satisfying that -> fi weakly and /W -> / in L 1 (D) for some f G B(D, [0, 1]). 
Then 

(i) the laws of (u^ N ') are tight in -D([0, oo), A4(D)) and limit points are continuous, 

(ii) for any 8 > 0, the laws of the occupation densities (u^ N '(5, t-\-8, x) : t > 0, x G D) are tight in C([0, oo), C(D)) 

where C{D) is the space of continuous functions with the topology of uniform convergence on compacts. 

Proof. The moment bound (|4.7[) allows us the uniform control on total mass moments: 

sup\u[ N \l)\A <C(p,T,/3, 7 )(l + (^ N \l)r). (4.11) 



E 



2G 



To control the moments of the approximate densities, the following simple upper bound for U W is useful. Triggering 
all the nutrient packages (uk '■ 1 < k < Kn) used in the construction of u^ N ' at time zero gives the process 



J2^=o u k,t- This satisfies 



,(N) 



(0, t, x) < u iN) (0, *, x) for all x G D and t > 0. 



(4.12) 



Moroever has the law Q D (2) + af(N) dx of a DW(D,7) process. 



Now fix <f> G Cq(D). Using the decomposition (|4.5p . we consider the increment (4>)—ui N ^ (<f>) \ for < s < t. 
This leads, via the total mass moment bounds and Burkholder's inequality, to the estimate 



E 



sup 

re[s,t] 



4 N \cf>)-uf >(</>) 



(N), 



<c(\t- s\ p/2 + N- id+1)p ^j for all N and s,t G [Q,T], 



(4.13) 



where C < oo depends only on P,j,4>,p,T and supjy/x^(l) (we suppress this dependence below). The term 
jy-(d+i)p comes f r0 m the size of the jumps in the Burkholder inequality, for a cadlag martingale M, in the form 
E[sup s<t \M s \p] < C p E[[M] p t /2 ] + C p E[snp s<t \D S M\ P ], The estimate (HT3)) gives, via a Chebychev inequality, 



sup 

se[t,t+N-( d + 1 i] 



< CN' dp/2 . 



By summing over the grid tj = jN ( <i + 1 ) and applying Borel Cantelli, we find that 

ui N) (<f>) - u[ N \(j>) < 2N~ 1/2 whenever < s, t < T and \s - t\ < N-( d+1 \ for all large N, P a.s. 



(4.14) 



We claim that 1|4.13[) is sufficient to imply that uS N ' (<ft) is tight in D([0, oo), R) and has continuous limit points. 



To see this we smooth the sample paths by setting u t 
u\ N \<j>)-4 N \cj>) <2N^\t- S \\\ 



(N) 



N 



"[t,t+iV-( d + 1 )]' 



oo sup 4 W) (1) for all N and s,t G [Q,T]. 



Then t — > w l '(4>) is continuous and 



r<T+l 



This together with P~T3|) implies that Sp^^) 



(</>) are tight in C([0, T],R). But (|4~T4f shows that the error sup t < T \u\ >{4>) 
(</>) has the the same continuous limit points as uf \<f>). 



(<t>)\ p ] < C\t - s\p/' 2 for s,t G [0,T]. Hence the laws of 

(0) | — > almost surely and 



7 W 



We now state a compact containment condition. Let D e = {x G D : d(x,D c ) > e}. Choose smooth 0£>.e : 
D — > [0, 1] satisfying (j>r>,e = on I?2e and = 1 on D\D e . We also suppose that is decreasing in e. We 
claim that for all S > there exists e(<$) > satisfying 



lim lim sup P 

5^-0 N-^oc 



sup U t 

t<T 



(4.15) 



We postpone the proof of this to the end of the section. This condition controls the amount of mass near the 
boundary and allows us to extend the tightness of u[ N \4>) above to all 4> G C®(D) by an approximation argument. 
It also is the compact containment hypothesis in [10] Theorem II.4.1, from which we obtain the tightness of 
in D([0,oo), A4(D)) and the continuity of the limit points. 

To establish the tightness of the occupation densities we use the following increment estimate. Under the 
hypotheses of Proposition [TU we claim there exist < ati < oo and finite C, depending on f3,-y,p,5,6',T,D and 
su Pjv fi {N) (l) + (f {N \ 1), so that 



E 



\u^(s, s ,x)-u^\6,t,y)\p <c{\t~ s r p + \ x - y m 



AN) 



(4.16) 
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for all 5 < s,t <T and x,y € D\ Dg'. This implies tightness of (t,x) — > (<5, t, x) in C([£, oo), C(D)), where 
C(D) is the space of continuous functions with the topology of uniform convergence on compacts. 

To establish (|4.16J) we use a Green's function representation for the density (note the density exists since it 
exists for each component part Ukt)- It is convenient to break the approximation into two parts 



,(N) 



u o,t + u\ where u 



(N) 



Uk,t-r k - 

keSt 



(4.17) 



For the process we may consider ^^(O, t, x) since the nutrient packages do not start at a singular initial 
condition. Indeed combining the Green's function representation (|3.19p for each U] i t _ Tk and summing over k S S t 
we obtain 

i m (0,t,x) = f f G^_ s] (z) {- 1 ui N \dz)+m^(dz,ds)] + Y, D ' u '( G ?o*-.)) (4.18) 

JO JD s < t 

where mS N ^ is the martingale measure associated to iS N \ We now consider increments of this representation. 
There are many terms, so we choose to illustrate the idea on one key noise term and on the jump term. We shall 
also state all the underlying estimates on the kernel Gp'^. First 

Gp%(*) < G [ll\( z ) < C \x ~ for t > and x,zeR 3 

and the analogous bounds G^£(z) < C(T)(1 V ln(l/|a; - z\) and G^Jz) < C{T) when t < T. Also 



for all t, t' and x.z <G D and 



\ G [o%( z ) - G"o%(*)\ < C(d, S,T)\x -y\(l + \x- z\'- d ) 



for t,t' < T and x,y G Dg, z S D satisfying \x — z\ > 2\x — y\. These estimates are explicit calculations when 
D = R d . One way to obtain them for general D is via a suitable coupling argument. 

Consider the term K(t, x) — J Q J G^ j '^_Jz)rh^ N \dz,ds). The time increment K(t + s,x) — K(t,x) can be 
split into two terms, defined by 



G^t+s-rp) ~ G» : l r] (z)mW(dz,dr), K 2 



t+s 



G^_ r] (z)m^>(dz,dr). 



Using a Burkholder inequality, and the kernel estimates above we find, when d = 3, 



E[\Kx \p] < C( P )E 



D 



u {N >(0,t, z) (\z - x\~ 2 A s 2 \z - x\- 6 ) dz 



p/2 



A comparison as in (|4.12p combined with the moments (|3.21j) shows that sup z E[uW{0,t,z)P/ 2 } < C{p,(3,T) when 
t < T and z £ D. This leads to E[\K X \ p ] < C(p, /3,T, D)s p / 4 . Similarly, again when d = 3 and t + s < T, apply 
Holder's inequality to find 



E[\K 2 \P] < C{p)E 
< C(p)E 



I(re [t,t + s])ui N \\x--\- 2 ) 
p/s 

I(re [t 7 t + s})ui N \dz) 



|x--r 8 / 3 ) 



3p/8 



< C(p)\s\ p/8 E 

< C( P ,(3,T,D)\ s \ p / s 



8Upfii W) (l)|*/ 8 
t<T 



u iN) {0,T,z)\x~ z\- 8 ^ 3 dz 



3p/8 
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In this second estimate we do not have the optimal power of s, but have shown how to rely only on the simple 
moments estimates for sup t<T u[ (1) and u^ N \0,T,z). 

The jump term in the representation (|4. 18|) can also be easily handled, since the sum of all jumps in the 
measure t — > ut is at most (3dx. For example, an increment in x is bounded using the kernel estimates above, when 
d = 3 and \x - y\ < 1/2, 



Y. IKu - [(:! ''j - : - D » u >( G \o'*-.]) 

s<t 

< p [ sup\G°'*(z)-G 

JD s<t L ' 1 



< C0 \x- z^ 1 + \y- z\- 1 dz + Cf3\x-y\ / 1 

J Bix^x-y^/ 2 ) J D\B{x,\x-y\ 1 / 2 ) 

< C(P,D)\x-y\. 



l dz 



Other terms in the Green's function representation of the increment u(Q, t, x), and in all dimensions d = 1, 2, 3, can 
be controlled by similar estimates. 

From (|4. 1T[) we have u^ N \S,t, x) = uo(5,t,x)+u^- N \5,t,x) and we have controlled the increments of u^ N \0,t,x), 
and hence of u^ N ^(S,t,x) = u^ N '(0,t,x) — u^ N \0,8,x). The increment estimates for u (S : t,x) are similar, except 
that one restricts to 6 > and uses (|3.20p . One also needs estimates on the fixed time measure u s to control the 
increments of uqj(G? q '®_ s \ For this, use the method of Lemma IIL3.6 in Perkins [TU] to show, for p G [0,2) and 
suitable Co = co(p) > 0, that E[exp(coUo,s{\x — < C(d, sup ra // n (l), S,p) < oo, which is sufficient. □ 

Proof of the compact containment ()4.15|) . We start by showing that 



sup £7 

N 



J D,e 



)M 



as e | 0. 



Indeed, arguing as in (|3.28|) we have 



Q"/' 1 [EW02V)A1] 10 ase|0. 



(4.19) 



(4.20) 



Moreover the continuity of the laws fx — > Q^} n on C([0, oo), M(R d )) ensures that the expectation in (|4.20|) is a 
continuous function of /i. Therefore the limit in (|4.20|) is uniform over the compact set JCd.l of measures supported 
in D and with total mass fj,(l) < L. By the coupling (|4.12|) and the simple comparison (I3.29|) we find 



E 



[U % T]{4>D,e) A 1] < Q%fi +fljmda [U [0 ,T](<l>D,e) A l] 



which leads to (|4.19j) . 

To control the supremum s\xp t<T u^ N \<f>D,e), we will establish a modulus of continuity that is uniform in N 
and e. Let G^4> denote the action of the Green's function as defined by Gf<f)(x) — f D G®' x (y)4>(y)dy. Then, 
arguing as above, we have 

sup sup fi (N) {G?<t> D ,e) -► as e -> 0. (4.21) 

N t<T 



For <f> £ B(D, [0, 1]), the measurable functions from D to [0, 1], define 



X t {N,<t>)=u\ N \cl>)-^ N \Gf(j>) and X t [N,4>)=N 



d+l 



We may extend the decomposition (|4. 5[) to time dependent test functions as in section 3.1.1. For smooth cf> 
compactly supported inside D, using the test function (s,x) — ► Gf_ s (j>(x) in (|4.5[) . over the interval s € [0,t], leads 
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to the Green's function representation for Xt(N, 4>). This representation in turn can be used to reach an increment 
estimate 

E[\X t {N,<t>) -X s (N,(j>)\ p } < C (\t - s\ p/2 + N- (d+1)p ^ for all N and s, t G [0,T], 

where C depends on p, T but not on <p or N . This requires (see Perkins |12j Corollary 5 for this argument) only a 
smoothing property of the Green's kernel, namely we use: 

\Gu' x (y) " G°' x {y)\dy < C{D)r l {u - t) for all < t < u. 



sup 

xeD Jd 

This estimate holds by direct calculation for the case where D = R d , and thereby when D is a box, since then 
the Green's kernel is a finite combination of reflected copies of the free space Green's kernel. For domains with a 
smooth boundary, it can be derived from Corollary 5 of Davies pQ . 

The smoothed version X t (N,<f>) then satisfies E[\X t (N,4>) - X s (N,<f>)\ p ] < C\t - s\ p / 2 for all N and s,t G 
[0,T]. The argument from Perkins [12] Theorem 1 then shows that for all 4> G B(D, [0, 1]), there exists a random 
c((f>, N, T) > so that 

\X t (N,(j>) -X s (N,(f>)\ < |t — s| 1/3 whenever s,t G [0,T] and \t-s\ < c{(j>,N,T) (4.22) 

and moreover that P[c(<p,N,T) < S) — > as S J. uniformly over N and <j> € B(D, [0, 1]). The point here is that 
the estimates on c(<p,N,T) depend only on total moment bounds and on ||</>||oo- 

It is straightforward to combine (|4.19[) . (|4.21[) and (|4.22[) to see that there exists e(S) > so that 



lim sup P 

S^o N 



supu[ N \(f>r> te u)) > 2S 

t<T 



0. (4.23) 



Indeed, consider the set 



supu 

t<T 



u ( t N) ((t>D,e)>28, c(0 Ae ,7V,T+l)>5 : 



If S > 1 and supjy sup i<r+1 f/,( N '(Gf ((/>D,e)) < f5 4 /8, then on this set we have sup t<x Xt(N, 4>D,e) > 5 and the 

modulus of continuity ensures that J Q X t (N, (j>D,e)dt > S 4 /A, and thence that J (<f>D,e)dt > S 4 /8. But 
this has small probability by (|4.19J) and Markov's inequality. So (|4.23p follows by first choosing 8 so that that 
P[c(4>d,(, N, T+l) < S 3 ] is small, uniformly in N and e, and then choosing e so that sup N sup t<T+1 fj,( N '(Gf \<j>D.e) < 
(5 4 /8 and sup^y E\u^^ ((f>D,e) A 1] is small. 

To remove the time smoothing in I4.23[) . we would like to apply (|4.14p . However, since <j>r> e does not lie 
in Cq(D), there is an extra exit measure term when deriving the increment estimate (|4.13|) . However the exit 
measure is non-negative and we can still deduce the one sided estimate u^(cf>) < ui N \(f>) + 2N~ 1 / 2 whenever 
< s < t < T and |t — s\ < N~^ d+1 \ for all large AT, P a.s, and this is sufficient to deduce from (|4.23p the desired 
compact containment estimate. □ 



4.3 Passage to the limit 



The aim now is to prove Theorem 1 1 5 1 below . showing that the approximations constructed in section |4T] converge 
in law to solutions of (jl.ip . The method is to pass to the limit in the approximate martingale problem (|4.5|) and 
confirm that any limit point of the sequence of approximations must be a solution to (|l.ip . Many steps in the 
argument are quite standard (see, for example, the arguments in Proposition II. 4. 2 in |10j). so we concentrate on 
those that concern the interaction term and the exit measures. 

Since we wish to let N — > oo, we again include this dependence in the notation. Assume the hypotheses of 
the tightness Proposition Q3J For any subsequence, we may choose a sub-subsequence where the approximations 
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and their occupation densities converge in distribution. By a Skorokhod embedding, we may choose versions of 
these approximations that converge almost surely. Without changing notation for the version and labelling the 
convergent sub-subsequence still as N, we may suppose there exists a limit process u, with continuous paths in 
A4(D) and a continuous field (u(s, t,x) : < s < t, x £ D) so that, almost surely, 



.(AO 



sup 

x£A,s£[5,t] 



sup 

s<t 

((5, s, x) — u(5, s, x) 



(0) - u s (4>) -> 0, for all 4> G C$(D) and t < oo, 



0, for all < 6 < t < oo and compact A C D. 



(4.24) 
(4.25) 



Note that u(s, t, x) must act as a density for U[ st ](dx) for < s < t. Since also s — > u(s, t, x) is non-increasing we 
may set v t (x) = f(x) exp(— u(0, t, x)) — lim s | f( x ) exp(— m(s, t, x)). 



Fix <fi G Cq(D). Consider the martingale problem (|4.5p for vS N ' at a fixed i > 0, noting that 



dD,(N) 



(0) = o. 



The terms u[ N \<p), fj,( N \(f)) and J* ui N \A(f> — ^<j))ds converge pathwise as N — > oo by (|4.24p . The error terms 
E^' N \(jj) and Ej. 2 ' N \(j}) converge to zero (in L 2 and pathwise respectively) using the estimates (I4.4[) and (|4.6[) and 
the moment bounds (|4. . To handle the key term J* ui (vi <f>)ds we approximate, using (|4.10[) . 



S s N \vi N U)ds = J «W (/We- u [».-]^V)d fl + £ t (3 ' JV) (^) 
where the error term is bounded by 



< ll^lloo / « W (0>M) sup<' w (a;)rfi 



(4.26) 



and m" t ' N (x) — /J exp(— uf^,] (iJ>x))dm]!. ,N (x). Doob's inequality, and the estimate (|4.9j) . shows that 



sup|m f (x)\ 



< 2N- 



E 



ds < N~ 



The comparison (|4~T2|) shows that E [(u {N) (0, t, x)) 2 ] < Q D (N) +0Indx [(U(0, t, x)) 2 ] . Applying Cauchy-Schwarz to 
P~26"]) . and using the the estimate ([3~22]) . we find that E(\Ef' N) (<j>)\) -> as iV -> oo. 

For 5 G (0, i) we write 



r/.s 



(/e^V^ 0) ds + E^ N \4>) + E^ N \4>) + E^ N \j>). 



ds 



(4.27) 



The error term Ef' N ^(<f)) is bounded by ||</>||oc / u< s N \^-)ds and converges to zero, uniformly in N, as <5 — > 0. The 
error term E^ ,N ^ (</>) is bounded by 

ft 



11011 



The approximate density Ur^i(^x) satisfies f D (u^Q^(ip x )) 2 dx < J D (u^ N \0, 5, x)) 2 dx. This, together with Cauchy- 



u{ N \u^.))ds = \\ 



AN) 



11 ^»(0,i, a :)«[3(^) ( i I 



D 



Schwarz as above, shows that E^ 6 ' N ' (<j)) converges to zero in L 1 , uniformly in N , as S — > 0. The error term E^' 1 ^ > ((f)) 



(4.JV), 
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is bounded by 

\\4\\°oJ* 4 N) (\f W - f\) = UWooJ « W (0, t, x)\fW(x) - f(x)\dx. 

Since are bounded we have — > / in L 2 (D) and combined with (|3.22j) we see that E[ i,N \<j)) — > in L 2 . 

Finally, the first term on the right hand side of (|4.27|) converges pathwise, for fixed S, to L u s (f exp(— u(S, s))<j>) 
as N — ► oo. This follows from (j4.24l4.25j) when / 6 C®(D), and for general measurable / we can approximate in 
L 2 by continuous / and control the error as for the error term (</>). The limit fji u s (f exp(— u(5, s))<p)ds is, 

by repeating the above approximations, close, for small 5, to J Q u s (v s (j>)ds. The conclusion is that the key term 
converges, in probability, to J Q u s (v s (f>)ds. 

We may now define a continuous process mt(4>), for <fi G Cq(D), by the formula (|2.1jl . The convergence of 
all other terms ensures that m[ N \<p) —> mt(<p) in probability. Moreover, standard arguments yield that m t ((j)) is 
a martingale with respect to cr{u s : s < t}, and with the correct quadratic variation (note that we have uniform 
moment control by (j4. 1 1 j) ) . Thus the limiting process (it, v) is a solution to the martingale problem for (jl.lj) started 
from (/x,/) with respect to its natural filtration a{u s : s < t}. We have almost proved the following convergence 
result. 

Theorem 15. Let D be a bounded domain. Suppose the approximation (u^ N \v^ N ^) has initial tig ' = and 
v o N ^ = — 1 satisfying that weakly and f^ N ^ f in L 1 (D) /or some / G £>(-D, [0, 1]). Then, 

(i) the laws converges in distribution on D([0,oc),M(D))xC([0,oc),M(dD)) to the limit Q®j >n ; 

(ii) for any T < oo, the law of the triple (u!^,, u^ ' , fj"^) on M(D) x M(dD) x L 1 (D) converges to the law 
of (C/[o,T], t/l 13 , /exp(-[/(0, T))) wufer 

Completion of the proof. The convergence of the exit measures can be deduced from the extended martin- 
gale problem for test functions G C 2 (_D). Indeed, return to the sub-subsequence studied before the statement of 
the theorem. Since (u, v) solves there is a continuous exit measure process u dD solving the extended martin- 
gale problem. Using an approximation argument by </>„ G Cq(D), it follows that m t (<f>) — > m t (4>) in probability 
for all <fi G C%(D). Choose <f> G C 2 (Z?) so that h = 4>\od is non-negative on dD. Then passing to the limit in 
the extended approximate martingale problem (|4. 5 j) for u'"', we have shown convergence of all but one of the 
terms, so that this last term u^ D '^ N '(h) must also converge in probability to uf D (h). The fact that t — > u[ N,> ' 9D (h) 
and t -► u<=> D (h) are non-decreasing and continuous imply, at least along a further subsequence and for all t, that 
sup s<t \ui N ^' dD (h) — u d g D (h)\ — > in probability. A further subsequence argument allows us to conclude the same 
for a countable dense set in C®(dD), and this implies that C([0, oo), M(dD)). The uniqueness 

in law of (u,u dD ) for solutions to (jl.ip implies the convergence of (u^ N ', u dD '^ N ') in part (i). 

Based on (|4.10|) . the proof of the L 1 convergence of part (ii) uses similar, but slightly simpler, tricks as above 
when dealing with the convergence of the key term J* ui N \vi N ^ (f>)ds, and the details are omitted. □ 

5 Proof of the decomposition results 

This section contains the proofs of the decomposition results stated in section |2"T21 In section HTT1 we show that the 
exit and occupation measures of the approximate discrete nutrient process, as constructed in section 14. 1 1 can be 
done in two stages. In section 15.21 five different examples of this two-stage construction lead to five decomposition 
results for the approximation processes. Passage to the limit is done in section f5. 31 
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5.1 A two-stage construction of the approximations 



The construction of the approximation process given in section 14.11 is a pathwise construction. The approximation 
is constructed as a deterministic procedure, which we call the basic construction, applied to fixed realizations of 
the non-interacting DW(D, 7) processes (u^ : k = 0, 1, . . .) and exponential variables (e^ : k = 1, 2, . . .). In this 
section we show that the occupation and exit measures from this construction can be built up in two stages, where 
each stage applies the basic construction to a certain set of variables. This is described in a somewhat abstract 
manner, and the reader might want to look ahead at one of the five examples in section \57Z\ to be convinced it is 
a natural idea. For the rest of this subsection we act pathwise, supposing a single realization has been fixed of the 
underlying variables. 

Integrating over t in (|4.1|) . or letting t — > 00, we obtain 

W[ ,oo) = w ,[o,oo) + u fe,[o,oo) and u d J? = + ^ uf^, (5.1) 

where S^ = lim^oo St is the set of labels of nutrient packages that are ever triggered. So the total occupation 
and exit measures are determined by the set S^ and the total occupation and exit measures of the non-interacting 
DW processes. Moreover for k £ S^ we have 

^,[0,00) 

while for the k 5^ the converse inequality holds. This is exactly the condition that 5^ is a fixed point of the 
mapping T defined in the abstract lemma below, once the choices A — {1, . . . , K} and 

e(k) = e k , f{k) = Uo,[o,oo)(Wj M(j,k) = Mj,[o,oo)(W for j,k e A 

have been made. 

Lemma 16. Suppose A is a finite set, and fix e, / : A — > [0, 00) and M : A x A — > [0, 00). Define, for B C A, 

T(B) = iaeA: f(a) + £ M (b, a) > e(o) I . (5.2) 

I b£B J 

Then there is a unique smallest fixed point S C A of T, that is T(S) = S and S is contained in any other fixed 
point. Moreover if B C S then T n (B) equals S for large n. 



Proof. Since M(a, b) > we see that if B C B' then T(B) C T(B'). So T™(0) increases to a limit S which must 
be a fixed point of T. Any other fixed point S' contains T™(0) for all n and hence contains S and the uniqueness of 
the smallest fixed point is clear. If B C S then T"(0) C T n (B) C T n (S) = S and the result follows since T n (0) = S 
for large n. □ 

It is straightforward to check that S^ is the smallest fixed point of T. Indeed, let k m be the label of the 
mth nutrient package to be triggered, and let r m be the time it is triggered; then the definition of r m shows, when 
r m < 00, that k m £ T({ki, . . . , fc TO _i}). Inductively, {fci, . . . , k m } C T m (0) must hold and letting m — > 00 shows 
that Soo is contained in the smallest fixed point of T. 

Example 0. To illustrate the use of this lemma, we give here the analogue of the comparison Lemma[S]for the 
approximation processes. Suppose that Ylk=i ^ k ~ f ana - ^2k=L+i ^ k ~ 9- Then running the basic construction 
on the processes (ilk '■ k = 0, 1, . . . , K) and exponential variables (e^ : k — 1, . . . , K) we obtain the approximation 
process (u, u dD , v) started at fi and with initial nutrient level / + g, and the triggered set S^ that is the smallest 
fixed point of a mapping T. We now set Uq = uq + J2k=L+i u k- Now running the basic construction on the 
realizations of (uq, u\,U2,..., ul) and the exponential variables (e^ : k — 1, . . . , L) we obtain an approximation 
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process (u, u dD ,v) started at fj, + Pgdx and with initial nutrient level /, and the triggered set Soo that is the 
smallest fixed point of a corresponding mapping T. Moreover, for A C {1, . . . , L}, 

T{A) = jfc G {1,...,L} : u ,[o,oo){i>k) + Y1 T U e A ) u j,[o,°°)$k) > efc 

fe G {1, . ..,L} : Mo,[o,oo)(^k) + ^V(j G {£ + 1, .. .,if})%[o, o)(Vife) > e fe 

= T(AU{L + 1,...^})\{I + 1,...,X}. 

Choosing A = and then iterating we find that T™(0) \ {i + 1, . . . , K} C f "(0) and hence, by Lemma [TBI that 
5oo \ {L + 1, . . . , X} C ^oo. Combining this with (|5.1[) shows that (tt[o j0 o), w^P) is smaller than (it[o,oo): "So 9 )- 

We now describe the two-stage procedure, which will build up the set iSoo of triggered nutrient packages in 
two steps. We suppose there is a splitting of each of the DW processes (uk ■ k = 0,...,K), and their exit measures, 
into two parts, and , each with continuous measure valued paths, and satisfying 

/ dD \ ( - . + dD- . dD, + \ ,r on 

(«fc,[0,oo),«fc,oo) = ( u fc,[0,oo) + U fc,[0,oo)' U fe,oo + M fc,oo ) ■ ( 5 - 3 ) 

In stage one we apply the basic construction using the processes (uj, : k = 0, 1, . . . , K) and the exponential variables 
(et- : k = 1, . . . , if) to create a process t>~) and a set of triggered nutrient packages. In particular 

>,oo) = >[0,oo) + U fc,[0,oo)- ( 5 - 4 ) 

We then define {to,t = u^ t + J2kes^ u k t an< ^' for fc G {1, . . . , K} \ S~ , 

Uk,t = u k ,t and e fe = e fe - WQ i[0iOo) (Vi fe ) - ^ wT [0)Oo) (^ fe ). (5.5) 

In stage two we run the basic construction on the processes uq an d ("fc : k G {1, . . . , iC} \ <S^), with the nutrient 
packages triggered using the values (ik ■ k G {1, . . . , K} \ S^). This leads to a second process (u + , v + ) and a 
second set of triggered nutrient packages <5>+ C {!,..., if} \ <S~ satisfying 



" " ' : " ' * W fc,[0,oo) + J2 U k,[0,co)- (5.6) 



Adding (|5.4p and (|5.6[) and comparing to (|5.1[) . we see that U[ 0jOO ) = itj n + , will hold provided that the 
equality 5^ = U 5+ holds. A similar argument shows that this equality is also sufficient to ensure that the 
total exit measures satisfy m^? = u^'~ + u^' + . 



To verify this equality we use the following second abstract lemma. With the choices 

e(k) = e k , f ± (k) = u± [0tOo) (4> k ), M ± (j,k) = «J [0iOo) (^fc) iorj,keA 



note that the definition of /, e, M in the lemma corresponds to the definition of Uq, efe in the stage two construction. 
The conclusion of the lemma then implies that the equality Soo = U <5>+ does indeed hold. 

For the statement of the lemma we denote the unique smallest fixed point in Lemma [TBI by S — S(A, e, /, M) 
to indicate its dependence. 

Lemma 17. Suppose A is a finite set, and fix e, f~ , f + : A —* [0, oo) and M~ , A/ + : A x A — > [0, oo). Set 
M = M~ + M + and f = f~ + f + ■ Then the fixed point S(A, e, /, M) can be broken into two subsets as follows. 
First let S~ — S(A, e, / _ , M — ). Then define, for a G A \ S~ , 

f(a) = }+(a) + M+(b,a), e(a) = e(o) - f~ (a) - ]T M-(b,a) 
bes- bes- 
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and write M for the restriction of M to (A\S )x(A\S ). Then e > and if S + is the fixed point S(A\S ,e,f,M) 
we have 

S(A,e,f,M) = S-US + . 

Proof. Let T (respectively T~ and T + ) be the map defined by (|5.2[) used for the definition of S (respectively S~ 
and S + ). The fact that e(a) > for a G A \ S~ follows from the fact that S~ is a fixed point of T~. 

Since f> f~ and M > M~ we see that T~(B) C T(S). (Note in particular 5" C T n (S~) for all n > 1.) 
Then (T") n (0) C T™(0) for all n and hence 5" C S. 

For B C A \ S~ we have 

T + (B) = aeA\S- : /(a) + ^ M(6, a) > e(a) 1 

I bG-B ) 

= \aeA\S- :/+(a)+ ^ M + (6,o) + ^ M(6,a) > e(a) - /"(a) - ^ Af-(6,o)l 

I beS- bGS bGS- J 

= JaG :/(a)+ M(b,a)>e(a)\ 

I bG-BUS- J 

= T(BUS-)\S~. 

Applying this first to B = and iterating, we obtain (T+) n (0) = T n (S~)\S-. The equality here uses S~ C 

Using the final statement of Lemma [TBI we may let n — * oo to obtain 5+ = 5 \ S~ as desired. □ 



5.2 Comparison theorems for the approximations 



Notation We denote the law of the approximation (ut,uf D : t > 0) constructed in section E3~T1 on the space 
D([0,oo),M{D)) x C([Q,oo),M(9£>)), by Q*f' 0n . 

We now give five examples of the two stage construction of the last section, each leading to a decomposition 
theorem for the discrete nutrient approximation process. In each example we describe the processes u^,u~^ used 
for the splitting ([S13]) . We let 

g- = a{u-} V a{S~} V a{u^I(k G S~) : 1 < k < K} 
be the information gained by observing the stage one process. 



Example 1. Decomposition for the initial condition /i. Suppose that \i = + Take independent 
DW(.D, 7) processes (uqA and (u^A, with initial conditions fi~ and fi + . Then uo t t = u^ t + UQ t is still a DW(D, 7) 
process. Using independent DW(D,j) processes (ut : 1 < k < K) started at Uk,o — ipk, we set 

ul t — Uk,t and ut t = for 1 < k < K and t > 

so that the splitting (|5.3|) holds trivially. The first stage leads to a process (u~,u , '" ") with law <9^- D /' 7 . The 
memorylessness property of exponentials implies that, conditional on <J - , the variables (e^ : k G {1, . . . , K} \ S^) 
defined in (|5.5p remain independent rate one exponential variables. The processes Uk for k G {0, 1, . . . , K} \ 5~ 
are, conditionally on (?~, independent DW(D, 7) processes. Note that — J2kgs~ ^ k * s &~ measurable. The 
second stage therefore produces a process (u + , u dD ' + ) which, conditional on Q~, has law Q N ' D '^ n , 
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Example 2. Decomposition for the initial condition /. We split the set of nutrient package labels into two 
parts {1, . . . , L} and {L + 1, . . . , K} and set /~ = J2k=i V'fci / + = Y^k=L+i V'fc- Then we form the splitting (|5.3|) 
out of the independent DW(D, 7) processes (u k : < k < K ) started at fi,ipi,,.., i\>k as follows: for t > 

_ _ ( u k . t for < k < L, + _ ( for < k < L, 

Uk ' t ~ \ ' for L + 1 < k < K, " fc >* ~ \ u k .t for L + 1 < k < K. 

The first stage of the construction produces a process (u~ ,u dD '~) with law Q 'f- ■ Note that — u = 

P^2k=L+i^kI{k £ S^,) and = X)fc£s~ V'fc are measurable with respect to the sigma field Q~ . Moreover, as 

above, the variables (efe, u k : k G {1, . . . , A'} \ S^) are, conditionally upon independent exponential variables 

and DW(D, 7) processes. The second stage therefore produces a process u aD ' + ) which, conditional on Q~ , has 

law Q N ff\ 
u o> v o 



Example 3. Spatial Markov property. We fix domains D C D + . We will choose the finite partition D + = 
UDj used in the construction of the approximation so that each Dj is either a subset of D~ or a subset of D + \ D~ . 

Then a subset, which we list as (ipk '■ 1 < k < L), will satisfy Ylk=i ^ k = fI(D~) an( i J2k=L+i = fI{D + \ D~). 
We apply the spatial Markov Lemma [8] to find processes (uk,u^,u^ : < k < K) so that {u k ■ < k < K) 
are independent DW(D + ,7) processes with initial conditions /x, ipi, . . . , ipK] : < k < K) are independent 
DW(D~, 7) processes with initial conditions [iI(D~), ip%,. . ., ipL,0, ■ ■ ■ , 0; and, conditional on u{u k : < k < L}, 

{ul : < k < K) are independent DW(L>+,7) processes with initial conditions Uqq — /iI(D^ \D ) + ^ tx^l-D+i 
\ D + for 1 < k < L and u^ = ipk for L + 1 < k < K. In addition, the Lemma [5] ensures the splitting 



+ dD- 



U k,0 = U k 



00 



(^fc,[0,oo),^ + ) = ("fc,[o,oo) + U fc,[0 lO o)> U £? ' \dD+ +U 9 k D ^' + ) 



for < k < K. 

This gives a suitable splitting as in (|5.3[) for the two part construction, where the exit measures are all evaluated as 
measures on dD + . The first stage of the construction produces a process (u~ ,u dD ,_ ) with law Qjjm-fjjm-)' 
The second stage produces a process (u + ,u dD+ ' + ) which, conditional on Q~ , has law Q N + D , where = 
Vl(D+ \ D-) + u^~f\ D+ and v+ = fI(D+ \ D~) + v~I(D-). 



Example 4. Decomposition for the birth rate /3. Fix (3 — (3~ + (3 + . Take independent DW(D,7) processes 
(u^ : < k < K) with initial conditions Uq~ = /i, Uq — 0, u^ — (3~tpk and u~£ — (3 + ipk- Set Uk = so 
that (uk '■ < k < K) satisfy the splitting (|5.3j) and are themselves DW(£>, 7) processes. Then the first stage of the 
construction produces a process u dD '~) with law Q^'f' 13 n ■ The second stage produces a process u dD,+ ) 
which, conditional on Q~ , has law Q N _i. D + 7 , where itj = /3 + (/ — w^) efe and = u^. 



Example 5. Decomposition for the death rate 7. Fix 7" = 7 + 7+. Lemma [TBI below ensures we can find 
processes (u k , w^T , : < k < K) satisfying the splitting (|5.3p and so that (uk : < k < K) are independent 
DW(D,7) processes with initial conditions n,ipi, ■ ■ ■ ,ipK', (u k : < k < K) are independent DW(D,7~) processes 
with initial conditions /i, ip\, . . . , ipx] and, conditional on er{u,7 : < fc < K} 7 (u^ : < k < K) are independent 
DW(D,7) processes with initial conditions (7 +u fc[ 0oo ) : < fc < K). Then the first stage of the construction 

produces a process (u~ ,u 9D '~) with law Q 1 ^'^' 13 ' 1 ■ The second stage produces a process (u + ,u aL>:+ ) which, 
conditional on Q~ , has law Q^ D ^ n , where = 7 + U[q ^ and = u^,. 

Lemma 18. Suppose j~ = 7 + 7 + . There exists a coupling of three processes: u a DWfj, D) process with initial 
condition fi; u~ a DW(^~ ,D) process with initial condition /i; and u + which, conditional ona{u~}, isaDW("f,D) 
process with initial condition 7 + "jq x y and moreover these processes satisfy the splitting 

(«[o,oo),«£P) = (>,oo) +«;,„).«»•■ +^' + ) ■ 
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Proof. Construct a coupling of two processes (u~,u) as follows. Let u~ be a DW(D,j~) process with initial 
condition /i. Conditionally on a{u~}, let u be a DW process on the space D x R, with spatial motion that is 
Brownian on the first component D and zero on the second component R, with annihilation rate 7 and with initial 
condition j + u r (dx)I(r > 0)dr. If we define 

ut(dx) =u t (dx x R) and uf A+ (,4) = uf (Dxfi) (A x R) for A C 3D 

then conditionally on er{u~} the process u + is a DW(.D, 7) process started at 7 + Ujq * and with exit measure 
u , ' + . Define measures / on [0, 00) x D and on [0, 00) x <9L> by 

/> CO /'CO /' 

I(fc< 2 )) = I-(hW) + I+(hW)= u-(hW)ds+ / ^(aOfiatda;, dr)ds, 

JO JO JDxR 

poo poo p 

I 9D (h^) = I aD '-(hW) + I dD > + (h&)= d s u 9D (h^) + / ^5 r ( JB )d.tt.(d ! r,dr) 

J JO Jd(DxR) 

for bounded measurable /i (2) : [0, 00) x D -> i? and h (3) : [0,oo) x dD ^ R. We claim that (I(dt,dx), I 9D {dt,dx)) 
has the same law as (Ut(dx)dt, dtU 9D (dx)) under Qf? n ■ This allows us to define the the required process (it, u 9 - ) 
using (J, and the required splitting follows from the definitions of / and u + . 

To establish the claim on the law of (I,I dD ) we fix smooth non- negative hs (x), hi 3 \x) that vanish for s > t. 
It is enough to check that E[exp(—I(h^) — I 9D (h^))] — exp(— fi{^ t )) where {<j> a : s € [0,i]) is the solution to the 
log-Laplace equation (|3.7[) . We start by calculating the conditional expectation E?[exp(— I(h^ 2 ') — I 9D (h^))\a{u~}]. 
Note the almost sure limits 

I+(h^)= lim V / /i2*(a:)fi.(da:,dr)ds= lim V / X K , s (hf )ds 

and 

Z 9A+ (M 3) ) = lim / / / r£[ ,,^ i) ^ 4 (a ; )d s « s (dx,dr)= lim £ / d s X|>f> ) 

where (X ky t ■ t > 0), for k = 0, 1, . . ., are defined by Xfc^A) = ut(^4 x ^r)) and are, conditionally on er{it~}, 
independent DW(D,7) processes with initial conditions Xk,o = l + u~ k h+x .. Using the Laplace functional (|3.6p of 

(X k ,X 9D ) we find 

E[e-' + ^-' aD ' + ( h(3 \{u-}] = lim exp(-7 + Ev ^Mt-*)) 



jV->oc 

e 



7 + Jiu-{<t>t-s)d8 



Then we complete the expectation as 

B [ e -J(h<'>)-/ BD (I.<">)] = J"q < (h< a >+7+*t_.)d.-/« d s «f"-(/ l < 3 )) ] 

which can be calculated again using the Laplace functional of u 9D '~). It is straightforward to check the 
required log-Laplace function is solved by <f> and this proves the claim. □ 

Remark. A more natural splitting for the final example is to let (uZ,ut) solve, for each k, the system 

dtu" = Au~ — 7~u~ + Vu~W~, 

dtu + = Au + + 7 + u~ — 71* + + vu+ VT + , 

with Uq — and with orthogonal martingale terms. Then uZ is a DW(/J, 7—) process and u k — uj, + is 
a DW(D, 7) process. The first stage of the construction produces a process (u~,v~) with law Q^'f 7 • The 
second stage produces a process which, conditional on Q~ , has the law of the approximation process 

with parameters /3,7, but with an extra immigration term 7 + u~. This leads to a corresponding decomposition for 
solutions with immigration to but we shall not make any use of it. 
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5.3 Completion of the proof of the comparison results 



This section contains the proofs of the decomposition lemmas SI [SJ [5] and using the comparison results for the 
approximations from section 15.21 with the convergence of the approximations established in section 15.11 

We start with a simple estimate controlling the death time r = inf{< : U t = 0} for the approximation processes 
on a bounded domain uniformly in TV. The change of measure argument given at the beginning of section \'d. 31 can 
be used to show the limiting death time has a tail Q^'f ' 7 [t > t] < C(/3,p, \D\, ^)t~ p for any p < 1. This argument 
does not apply directly to the approximations since they are not absolutely continuous with respect to a DW 
process. However, the argument below follows the same main idea. It is combined with a first moment argument 
that leads to a suboptimal bound, but which is sufficient for our needs. 

Lemma 19. There exists C{(3) < oo so that for all t > e and N > 1 

Q iW,7 [T > t] < cw ( (/i (i))i/2 + (/) !)) ln(r 1. 



Proof. Take an approximation process u with law Q^'f* ■ Taking (f) = 1 in the martingale problem (|4. 5|) 
and X s = — 2/s, so that X s = — we find for fixed t > and s e [0, t) 

de -x t _ s u s (i) = e -A t _.u.(i) L (i)(A t _ s + 1a?_ s - l\t- s )ds + \ t - s du d s D (l) - (3A t - s U s {v s )ds 



^L' ;j (e _At - stls(1) ) + A t _ s e _At - sUs - (1) £) s it s (l) + martingale increments 



s<t 



> -f3X t - s e 



u s (v s )ds + martingale increments. 



Let m(dx,ds) be the continuous martingale measure constructed from the martingales m t (0) in (|4. 5|) and, as in 
the change of measure arguments, set M t — J Q J v s (x)m(dx,ds). Then Ito's lemma, using the cross variation, 

d k(-/3M),e- Af — Ml = /3\ t _ s e- x *-* u *W£ s (-{3M)d[M s> m s (l)} s 

J S 

= f3\ t - s e- x *-* u «W£ s (-(3M)u s (v s )ds 

shows that s — > £ s (— /3M)(1 — exp(— A t _ s it s (l)) is a non-negative supermartingale on s £ [0, i). Taking expectations 
at s n 1 1, we find 

E [£ t (-(3M)I(T > t)] < lim E [f s „(-/3M)(l - e - A *-=™ u *™M) 
< (l-e-^A < 2fj,(l)/t. 
Let gk = inf{s : u s (v s )ds > K}. Then, using Cauchy-Schwarz, and v s < / < 1, 

(P [t > t, ctk > t]f < E [I(t > t)£ t (-0M)\ E [l(a K > t)£^{-f3M)} 

~I{a K > t)£ t {+l3M) e P 2 ti u °^ ds 



t 

To control o~k we use a crude first moment bound. From (14.81) we obtain 



~ t 



(5.7) 



£[(u t , 1)] = (/, 1) - / / E[v s (x)uM]dxds = (/, 1) - / P[ Us (^ s )]d S 

JO JD JO 

(using the fact that v is constant on the partition sets Dj). Then P\ok < t] < (f,l)/K by Markov's inequality. 
Combining this with (|5.T[) and the choice if = clni, for small c = c(/3), leads to the desired bound. □ 
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We may characterize the total occupation and exit measures via the Laplace functional, defined by 

for v x <EM{D),v 2 £ M{dD), < G C* h °(L>) and < V G C%{dD). Fix /i G M(-D) and / G ^(.D), for a bounded 
domain _D. Choose approximations processes with initial nutrient levels 1 > f( N ' — > / in L 1 and weakly. 
Theorem [T5] (ii) implies that 



lim Qj#5& [$0^(C/[o,oo),^ D )] = Hm lim [*^[o,t], ^)] 

iv— >oo " ,J iv — *00 i — ►OQ ^ iJ 



4 D ) 

T dD\ 



= T lim Jirn^ Qj#™ > [^(^[o,T], C/| D )] 



(5.8) 



The interchange of limits is justified by the above lemma, which gives uniform (in N) control on Q N ^' [Ut = 0], 



and the fact that U\ 



[o, 



U [0 ,r\ and [/£f = U%P on the set {U t = 0}. 



Proof of Lemma[6l Returning to Remark in section ISTTl we see that, taking approximations f^ N \g^ < 1 
converging to /, g as above, and noting that g^ N \x)dx — * g(x)dx in A4(D), 

Passing to the limit as above we find that 

q d J? 9 [nAUio^^^)] > q^% xJ [^(^.oo),^)] ■ 

This inequality between the Laplace functionals implies the desired stochastic domination. 

Proof of Lemma [4] part (i). Consider the decomposition u~ , u + given in section [5T2l example 1, with initial 
conditions jj, = [i~ + jjT and — > / in L 1 (D). Then pathwise 



an, 



£f 7 [^(Efaoo),^) 



lim E 

N-*oo 



= lim lim 

= lim Q A/3 ' 7 



As above, the expectation of the left hand side converges, as N — > oo, to Q 
law of (u + ,u + ) conditional on the limit of the expectation of the right hand side can be written as 

^<t>,ip(U[o,T],U^ D ) Q^fj y e -u(,o,T) [®<f>,ip(U[ 0tOO <),U^)] 



(5.9) 
Using the 



The interchange in limits follows as above, using = on {w^ = 0}. The N —> oo limit holds using the 
convergence of {u7 Q r ,, u dD ~, ) given by Theorem [T5l (ii) and (|5.8[) above. The final identity between Laplace 
functionals implies the desired result. 

Proof of Lemma [5j The proofs of both parts of this lemma follow closely the argument used for the proof 
of Lemma [4] part (i) above. For example, for part (i) take the decomposition u~,u + given in section T5.2I example 
4, so that (|5.9j) holds, and with initial conditions /i and — > / in L 1 . Then argue as above, noting that 



lim E 

N^oo 



[0,T]' U T 



~)Q 



0+U-v-) 



$<£,V>(^[Q,T]> ^T°) ^ i g+/(l- e -Pto.T))dx,/e- p C(i,T) (^[0,oo) > k^)] 
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using the fact that if v T converges in L l (D) then (3 + (f — v T )dx converges in A4(D). Part (ii) is entirely similar- 
starting with example 5 in section ["5.21 



Proof of Lemma We can follow the previous argument closely, but we write this out since we have 
to be a little careful about the two domains. To indicate that we use the joint Laplace functional acting on 
measures on D + and dD + , we denote it by Take the decomposition u~,u + given in section [5721 example 3, 



that 4J^(U[o,«.),«Sf ) = $ ^(>,oo) +>,oc)> 



,dD- 



,ao 4 



converges, as N -> oo, to Q^j' '' 1 ^^ U l '°°)' U ' 



rdD 



\dD+ + u ScT ,n ~)- Taking expectations, the left hand side 

For this proof only, set fi~ = fiI(D~), fi + = fiI(D + \D~), 

f- = fI(D~), /+ = fI(D+\D-) and fW ~ = f( N )l(D~), = f^I(D+ \D~) The limit of the right hand 

side can be evaluated as 



lim E 

N^oo 



,dD-,-\ \ n N.D + .,0, 1 
l oo \dD+ iV , an-,-, „ fAn , . 



lim lim E 

iV->oo T^oo 



lim D ,Pn 



^0>O U [O,T]' "T 

*U(U[o,T h W 



8£> - ,- 



3D" I 



' At + +^ D_ l D +,/+ + /e- c7 (°' T »/(£'-) 



)-D + ,/3,7 

'M++17S, 13- ,/++/e- c/ (°.°°)/(_D-) 



$+(c/ [0 , oo) ,c/£ D+ ) j 

*^(^[0,oo),C/^ + ) 

*£*(tf[o,oo),E& D+ ) 

+ ) 



The limits above follow as in the previous examples, once we have established that u^ D ~\qd+ converges in 
M.{dD + ) and u^P '~\d+ converges in M.(D + ). Both of these follow from the convergence of u^P '~ in Ai(dD~) 
under the hypothesis (|2.6[) . since this hypothesis ensures that the limit law of U® D under Q D _ '^l 1 does not charge 
the discontinuity set S = (dD~ n dD+) n (dD~ \ dD+) (check the first moment of U%P~ (S) is zero). 



Proof of Lemma [4] part (ii). To allow this case to follow from the same argument as the earlier examples 
we need a slight improvement in Theorem [15] part (ii). Taking the two stage construction from section [5T21 example 
3, we need to consider two parts of the approximate density, namely 

W = ^2 TpkI{k<£St) and u t + = ^ ip k I(k£Sf). 



Kk<L 



L<k<K 



Now consider a sequence of models indexed by N in which — > / ± in L X {D) and where / = / + / 



+ 



Then Theorem [151 part (ii) can be extended to show that the law of the quadruple (u} Q j,yU T 



(N) dD,(N) _(/V),- ~(iV),+ 



, Urp , Urp 



on M(D) x M(dD) x L^D) x L^D) converges to the law of (J7 [0)T] , t/| D , /" exp(-f7(0, T)), /+ exp(-C/(0, T))) 
under Q^f n - This can be shown as in Theorem [15] once one notes that the analogue of (|4.8p holds for v^, namely 
that there are martingales ' such that 



vf(x) = f(x) 



vf(x)u s {ip x )ds + m v t '^{x) 



v,±. 



The initial conditions for stage two of the construction can be expressed in terms of v as Uq = /3(f + — v^)dx 
and Vq = + and the passage to the limit then follows the same lines as previous examples. 



6 Life 



6.1 Embedded oriented percolation processes 



In the proofs of possible life, the main step is to construct a discrete one-dimensional oriented percolation (OP) 
(ui(j, k) : (j, k) £ C) where C — {(j, k) : j,k £ Z, j > 1, j + k is even}, so that life occurs if the percolation is 
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supercritical. The variables ui(j, fc) will be denned in terms of the exit measures for solutions to on a series 
of boxes. We will choose a parameter L > controlling the length scale. In the grid of points x^ k — (3jL, 2kL) in 
d = 2 or x^ k — (3jL,2fcL,0) in d = 3, the xi direction will play the usual role of time, and the xi direction the 
role of space for the comparison OP process. For this argument we therefore need to be in dimension d > 2. 

We will freeze the mass of the solution as it exits the boxes D(n) = (SnL, 3nL) d for n = 1, 2, . . . We define 
exit measures it^ 9 inductively. Choose an initial condition fj, supported in D(l) and let (u 1 ^ 1 ) be a solution 
on -D(l) started from (/x, 1). Suppose that {u J ,v 3 ) have been defined for j = 1, ...,n. Then, conditional on 
o\u? : j = 1, . . . ,n} we let (ii" +1 , t>™ +1 ) be a solution on D(n + 1) started from 

..n+l _ n,dD{n) ari A f T \ _ J W £> 0*0 X € D(n), 

u -Uoo , ana w j a, € D(n + l)\D(n). 

Repeated use of the spatial Markov property in Lemma[7]shows that u^ D has the same law as the exit measure 
U^ D(n) under the law Q^ ] ' 0n . 

In the hyperplane x\ — we define the box II — {x : |x| < L, xi = 0}. Define, for j, k G C, 

, m J 1 if^(^ + I,)>M, 
y otherwise. 

Set w(0, k) = if fc ^ and w(0, 0) = /(/i(Ji) > M). Define for (j, fc) e £, 



fc) 



1 if a)(j - l,fc - 1) = Q{j - l,fc + 1) 

Q(J, fc) otherwise. 



We aim to show, for suitable choice of the length and mass scales L and M, that (u>(j, fc) : (j, fc) G £) is a 3- 
dependent oriented site percolation with density at least 1 — e. We recall the definition of this from Durrett [1]: 
whenever (j n , fc„) G £, for 1 < n < iV, satisfy |j„ — Jm| > 3 or |fc„ — k m \ > 3 for all n 7^ m, then 

P[u(j n , k n ) = for 1 < n < N] < e N . (6.1) 

We write (0, 0) — > (j, fc) when there exist a sequence of points = fco, fci, . . . , fcj-i, kj = fc so that so that \k m — 
fcm-i| = 1 for < m < j and satisfying u;(m,k m ) = 1 for 1 < m < j. Then let Co = {(j, k) : (0,0) — > (j, fc)} be 
the cluster of sites connected to the origin. This definition is sufficient for the key property of percolation to hold: 
Theorem 4.1 of [4j shows that when e < eo = 6~ 196 then P[|Co| < 00] < 1/20, where |Co| is the cardinality of Co- In 
particular, if the initial condition is such that w(0, 0) = 1, that is > M, then 

Q^tpBDW 0] = p[^(n) ^ Q] > p[ | Cq | = ^ 

is bounded away from zero, and by Lemmas 1111 and 1121 possible life occurs. 

It remains to check the hypothesis (|6.1[) of being a 3-dependent oriented site percolation. We do this inductively. 
By conditioning on Qm = o\v? : j = 1, . . . , M} (and setting Qo to be the trivial sigma field) it suffices to show an 
estimate of the form 

P[u(M + 1, k n ) = for 1 < n < N | Qm] < e No . 

whenever |fc„ — k m \ > 3 for all n 7^ m. If ui(M,k n — 1) = oi)(M,k n + 1) = for some n then the conditional 
expectation is zero. So we may also restrict to the set where, for each n, there exists k n G {fc n — l,fc n + 1} for 
which ui(M, k n ) = 1. We let D(M, n) be the box x L M ~ + D 3L . Note that the boxes D(M, n) for n = 1, . . . , N are 
disjoint and contained in D(M + 1). 

We now apply the spatial Markov property Lemma[7]to the pairs of domains D~ — D(M, n) C D(M +1) = D + , 
for n = 1, . . . , ATq in succession. Since the domains D{M, n) are disjoint, this amounts to constructing, conditional 
on aju- 7 : j = 1, . . . , M}, independent solutions (u M+1 ' n , v M+1 ' n ) on D(M, n) with initial condition 

M+l,n — M,dD{M) T(r)(]\/[ and x _ f ^(x) X G D(M, n) fl D(M), 

Uq l{U{M,n)), and t, j x G D(M, n) \ D(M), 
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for n — 1,...,Nq. Then an exit measure u^ +1 ' dD ^ M+1 ^ can be constructed by running one further solution, 
conditionally on these No processes, that starts with the combined exit measures of (u M+1 ' n : n = 1,...,Nq) 
on their domains and any part of u^' 9D ^ that is unused, and run until the exit from D{M + 1). This leads 
to a construction of an exit measure u^, +1 ' dD ^ M+1 ' that, by the spatial Markov property, is equal in law to 
Um +1 ' dD ^ M+1 j but that is pathwise larger than the sum of the exit measures X)^=i u^ +1,n,dD restricted to 
dD(M + 1). This in turn implies, since the solutions (it M+1,n , v +1 ' n ) are conditionally independent, that 

P[lo(M + 1, k n ) = for 1 < n < N | Gm] < J] P[^ +1 ^ 9I " M '"»(^ (+1 ,„ + I L ) < M \ Q M ] 

n=l 

By translational invariance, and the monotonicity in the initial measure, the condition (|6.1[) will then follow from 
an estimate of the form 

Q D ^[U^{xl ±l +I L )<M]<e 

whenever / > I(x\ > 0) and fj, is supported on 1^ with ^i(l) > M . (6-2) 

By monotonicity in / it is enough to consider just the case / = I(x\ > 0). We shall check this block estimate in 
the next two subsections for two different regions of parameter values 7). 



6.2 Proof that liminf^o ^{(3)/(3 2 > in dimension d = 3 

After applying the scaling in Lemma[3]with the choices a = c = L/3 _1 ,6 = L 2 f3~ 2 ,e — 1, solutions to with 
the parameter values (/3,L _2 /3 2 ) in d = 3 become solutions to the rescaled equation (|6.3[) studied in the following 
lemma. Undoing this scaling, the lemma shows that the block estimate (|6.2[) holds for solutions to when 
fi G (0, 1] and with the choices 

7 = L,7 2 (e)/3 2 , L = L o {e)0-\ and M = ^(e)/3" 2 . 

Choosing e = 6$, the value from in section T6.ll that ensures the OP is supercritical, the OP comparison implies 
possible life for the parameter values (/3, 1/^~ 2 (eo)/3 2 ), and thus liminf ^{[3)/(3 2 > L^~ 2 (eo) as (3 — * 0. 

Lemma 20. Consider solutions, in d ~ 3, on the domain D3 = (—3, 3) 3 to the equation 

{ d t u = Au + L 2 /3~ 1 uv — u + ^/uW, 

\ d t v = -Lf3- l uv, v (x) =I(xeH), { ' ' 

where H = {x : x\ > 0}. For any e > there exists Lq = Lo(e) > 1 so that, whenever (u, v) is a solution to i6.3\) 
with parameters Lq and (3 S (0, 1], and with an initial condition fi supported in I\ = {x : x\ = 0, \x\ < 1} and 
satisfying > 1, 

P[u^'(xi,±i+h)<l] <e. 
where x\ t \ = (3, 2, 0) and X—\,\ = (—3, 2, 0). 



Remark. The intuition for the lemma comes from taking L large, ignoring for a moment the diffusion, killing 
and noise terms and reducing to the ODE system dtx — L^fi^xy, dty — —Lj3~ 1 xy, where, for L large, the solution 
converges fast to = xo + Lyo. Thus the nutrient leads to a large build-up of mass, which in turn leads to a high 
probability that the exit measure is large. The corresponding scaling in dimension d — 2 fails to produce this large 
build up of mass. 

Proof. We first construct a DW process u~ on D 3 with law Q? 3 ' 1 . Then, conditional on a{u~}, we construct a 
solution (u + ,v + ) to (|6.3p on D3 with initial conditions 

u+ = LI(x e H){l-e- L ^ lu '^°°^)dx and v+ = I(x € H ) e - L ^ lu ' {0 ' QO ' X \ 
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Lemma after a suitable linear scaling, shows that 



,dD 3 



/£? 3,+ has the same law as the total exit measure 



it^P 3 of a solution to (|6.3|) started at /i. Lemma [5] also implies that the exit measure u^ 3 ' + is stochastically larger 
than the exit measure U^ 3 where U has the law Q D + :1 of a DW process, conditionally on a{u— }. Choose smooth 
h : <9Z?3 — * [0, 1] so that {h > 0} = xi,i + I\ (the case of x\ t -\ is symmetric). Then 



E 



< E 




< E 









Q 



D 3 ,lr 



Hh)] 



= E 



where, using the DW Laplace functional from Q3.8|) . Aw = |ur + w on £) 3 and w = h on dD 3 . Let Co = inf{w;(x) : 
x G Z?2 H if}. Then a comparison argument shows that Co € (0, oo). Using the initial condition of Uq we continue, 
for j3 < 1 and L > 1, 



E 



-<? 3 (zi,i+/i) 



< Qu 3 ' 1 



cxp ( —L I [ 1 — e 

ffn.D s 



-Lp- 1 U(0,oo,x) 



w(x)dx 



< 



exp [ -Lcq I 1 1 — e 



-£/(0,oo,*)^ rfa . 



p(L,n). 



We will show that p(L, /.t) converges to zero, as L — > oo, uniformly over as stated in the lemma. Then a Chebychev 
argument finishes the proof. 

It is not difficult to see that if p(l) ^ and \x is supported on I\ then U[ ,i](i>2 n if) > 0, Q^ 3 ' 1 almost 
surely. It is possible to argue this directly from the Laplace functional, but it follows also from results in the 
literature. Indeed the projection of the process onto {22 = • ■ • = %d = 0} is a DW process started at cS , so 
one can use absolute continuity and the strong law of the logarithm for the closed support shown in Tribe (|14j). 
Since x — > {7(0, 00, x) is lower semi-continuous, the set {x G D2 H if : U(Q, 00, x) > 0} is open and almost surely 
non-empty by the above. This implies that p(L, fi) — > as L — > 00 for each non-zero p. Lemma 0] shows that 
fi — > p(i, /i) is non-increasing, so we may restrict to the set M.\(I) of measures supported in I and with total mass 
one. Note that Mi(I) is compact in the weak topology. Since L p{L,p) is also non-increasing, the required 
uniform convergence follows by a Dini argument if we can show that p — > p(L, p) is upper semi-continuous. This 
will be true since p(L, p) is the decreasing limit of p(L, /z, e) as e J, 0, where 



p(L,p,e)=Q° 1 



cxp —Lcq 



D 2 nH 



1 



D -C/(e,e" 



fix 



Increment estimates (similar to those in (|4.16p show that the laws Q^ 1 (U(e,e 1 ) G df) on C(D2,R) are tight, 
uniformly over /1 G A^i(i). This in turn implies that the map [i — > Q I j ) ll {U{e : €^ 1 ) G fi/) and hence also — > 
p(i,/i, e) is continuous, which finishes the proof. □ 



6.3 Proof that lim^^ P~ x ^{0) = 1. 

We fix k G (0, 1) and show there exists /3o = (3q{k) so that possible life occurs whenever f3 > /3o and 7 = k/3. Since 
^(/3) < /3 this implies the desired limit. 

Lemma 21. Fix n G (0, 1) and smooth 77 : D3 — > [0, 1] compactly supported inside D3 n {x : xi > 0} and so that 
77 = 1 on D2 ("1 {x : Xi > 1}. Consider, in dimensions d>2, solutions on the domain D3 = (— 3, 3) d to the equation 
for a DW process given by 

d t u = Au + 9 i^j «+ VoUW. (6.4) 

Given e > we may choose 9o(e),o~o(e) > so that for any solution u to J6'.^[ ) with parameters fJo and 6$, and 
initial condition /1 G Mi(h), 

> I - e 



uf D *(x lt ±i+h) > 1 



where xi.±i = (3, ±2, . . .) . 
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Proof. Choose h : dDj, — ► [0, 1] so that {h > 0} = xi,i + h (the case xi,_i being symmetric). Using the Laplace 
functional from (|3.6j) and scaling, we have E[exp(—uf Ds (h))] = exp(— /i(<^>i)) where 



&0 = A0 + 0(±±£t7-k)0-§^ on£> 3 , 



§0 2 

h on 9Z? 3 for t > 0, and </>o 



on D 3 . 



(6.5) 



As cr — * the solution <f>\ converges to 4>i where <f> is the solution of (|6.5p with a = 0. Since </> is decreasing in the 
parameter <r, the convergence is uniform over x G Ji. 

We claim that ^ (a;) — > oo as — > oo and that the convergence is uniform over x € I\. One way to see this is 
via the Feynman-Kac representation for <f>, namely 

fa(x) = E x \h(X T )I(r < lJe'^C^"^)-*)* 

where X is a Brownian motion (with generator A) and r = inf{i : X t £ dD^}. Note that if X s S -D2 H : Xi > 1} 
then (^^T](X S ) — k) > (1 — k)/2. Now we use the support theorem for Wiener measure. There is an open set of 
paths O x C C([0, 1], where / € O x satisfy fo = x, f exits D 3 before time 1 and Jq (^j 1 r](f s ) — K)ds > (1 — k)/4. 
Then ^(x) > exp(0(l - k)/4)E x [/i(X t )J(X € O x )]. Moroever O x can be chosen so that E x [h(X T )I(X. € X )] is 
bounded below uniformly over x G /1, establishing the claim. Then we may choose #o( e ) an d then <To(e) so that 



.0D 3 



(xi 



A) < 1 



< P 



< e e 



u? D3 (/i) < 1 



-M(0i) 



by Markov's inequality 



< ee + - by the suitable choice of a (e) 

< e by the suitable choice of {e) 

which completes the proof. 

Lemma 22. Consider, in dimensions d = 2 or 3, solutions on the domain D3 to the equation 

d t u — Au + 8u(v — k) + \J~ou W, 
dtv = —duv, vq = I(xi > 0). 



□ 



(6.6) 



Given e > we may choose 0i(e), o~i(e), <5i(e) > so that for any solution {u,v) to \6. 6\) , with parameter values 9\, 
o~\ and 8 < 81, and with initial condition \x supported on I± satisfying /x(l) > 1, 

P[u d °°(xx t±1 +h) > 1] > 1-e 



where ix,±i = (3> ±2, . . .) . 



Proof. By the monotonicity in /1 given in Lemma IH it is enough to consider /i € .Mi(ii). Fix #1 = 9o(e/2) and 
°"i = (J o(e/2) from the previous lemma. Let 77 be as in the previous lemma and let K be the closed support of n. 
We consider a solution to (|6 . 6|) on path space and set 

t§ = inf{t : 5C7(0,t,x) > (1 - «)/2 for some x € 

Note that V(t, x) = exp(-SU(0, t, x)) > 1 - (1 - k)/2 for i < tj and x G If. This implies that 

1 -|- K 

V(t, x) — n > — - — rj(x) — k for t < Tg and x G D§. 

Informally upto the time Tg the process dominates a solution to (|6.4p . More carefully, we may expand the process 
cxp(— u s (cf)i_ s )— u® D3 ((j) s )) over the interval s € [0, r^Al], with as in the previous lemma, to obtain an upper bound 
on the the Laplace functional ^[expt-uf 153 (h))I(r s > 1)]. This shows that P[u^ 3 (xi )± i + Ji) < l,r 4 > 1] < e/2 
for all fiE M.\{I\) (and for any choice of 8 > 0). 



44 



It remains to estimate P[t$ < 1]. Informally, we can control this by choosing S small and the fact that the 
process is dominated by a DW(£>3, —6\) process. Rather than use the natural pathwise comparison, which is messy 
(although possible) to establish in this non-Lipschitz setting, we exploit a change of measure argument. Let be 
the law of (u,u 8D3 ) on path space flr> 3 of the DW process satisfying dtu — Au + y/ayuW starting at /i. Set 



M t = a^ 1 9 1 J^ J (e su( -°' s ' x ^ - k) Z{dx,ds) 



and let £t(M) be the associated stochastic exponential. Arguing as in section [372! under the the measure dP\u t 
£t(M)Q^, the process U is a solution to (|6.6[) . Then 



P[rs < 1] = Q„ [£ ti ai(M)I(t s < 1)] < (q„ [^ A1 (pM)e ((p - 1)/2)[M] ^ A1 ]) 1/P (Q^r s < l]) 1,q 



(6.7) 



for a dual pair p^ 1 + q~ x = 1. Note that [M] t < (j^ 6\U\qa{1) for t < t§. Reversing the change of measure, under 
the law £(pM)dQ ll the process has a mass creation at most at rate p9\ up to time t$- For small enough 6 > the 
exponential moments Q® 3, ~ p61 [exp(0£/[o,i] (1))] are finite, and bounded uniformly over /i € A1i(/i). Thus we may 
choose p = p(e) > 1 so that the first term on the right hand side of (|6 . T|) is bounded uniformly over fj, € A'fi(Ii), 
and it remains to estimate Q^Tg < 1]. Note that 



\n < 1] < 



min{l, 25(1 - k) -1 sup 17(0, 1, x)} 

xeK 



:=p(S,im). 



The continuity of the occupation density U(0, 1, x) for x £ K (the finite spead of support ensures that the solution 
does not hit K immediately) implies that sup x&K {7(0, 1, x) is finite, and so that p(S, (i) J. as 8 J. 0. But p(^, <5) is 
the deceasing limit of p(5, /i,e) as e J, where 



p(5,n,e) = Q^. 



min{l, 25(1 - kT 1 sup [7(e, 1, a;)} 



Arguing as in the proof of Lemma 1201 the functions p(S,fi,e) are continuous in fi and hence p(S,fi) is upper 
semicontinuous, ensuring that the convergence p(S,fi) — > is uniform over Mi(Ii). This allows us to choose 8% to 
ensure that P[ts < 1] < e/2 for all /i <E Mi(Ii) and completes the proof. □ 



Now chooose e = eo the value that ensure the OP process is supercritical. Using the scaling lemma |3] for a 
solution to (jl.ip with the choices b = #i(eo)/3 _1 , a ~ b 2 , c d — abai(eo) and e = 1, the scaled solution u solves (|6.6[) 
with the parameter values 9\, o\ and 5 = 9\(3~ 1 . Choosing (3 sufficiently large that 5 < 5i(eo) we see that the block 
estimate (|6.2p holds, ensuring possible life. 

Remark. The correct large (3 asymptotics can be seen from the following formal argument. Choose 7 = [3 — 
p2/(6-d) for some q > q Applying the scaling in Lemma [3] with the choices a = /H^" 2 )^ 6 ^ , b = p- 2 / {& - d \ 
c = (3~ 1 /( 6 ~ d, > and e = 1 we find the scaled equation solves 

f d t u = Au + / 9(*-«0/(8-«0« (v - (1 - 6»/3-( 4 - d )/( 6 ~ d ))) + V^LU, 
\ Sju = -p-^-^/^-Vuv, v = 1. 

Letting v = 1 — /3 _ ( 4_d )/( 6 ~ d ){i we may rewrite this equation as <9tw = Au + 9u — uv + \fu~W and d%v = u(l — 
p-(4-d)/{6-d)fiy This suggests that as — > 00 the system converges to the one parameter system 

J <9 t u = Au + 9u — uv + y/uW, {f\S.\ 
\ d t v — u, v = 0. 

There is monotonicity of total occupation and exit measures in 9 and one expects that there is a critical 9 C so 
that solutions to (|6.8p die for 9 < 9 C and may live for 9 > 9 C . It is therefore reasonable to conjecture that 
hm^—Kx, f3~ 2 /( 6 ~ d \f3 — ^(/3)) = 9 C . There is no obvious line of proof for this conjecture, since it is not obvious 
that a life or death block estimate will hold for 9 close to 9 C . However such block estimates should hold for 
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sufficiently large and small 9 leading to upper and lower asymptotics of the same order. However note that the 
block construction we have presented in section T6.2I uses initial nutrient condition / = I(xi > 0), and the above 
intuition does not apply as such. One needs a more careful block construction where the nutrient level is controlled, 
and this in turn seems to require a stronger version of the spatial Markov property. We will present these details 
in a subsequent paper. 

We remark that a similar one parameter scaling limit should hold for the small (3 asymptotic. Indeed letting 
/3 — ► in (|6.3|) suggest that the small j3 behaviour when 7 = L~ 2 (3 2 should approximate the one parameter system 

d t u = Au + L5 a M — u + ^fuW where <j{x) = inf{7 : {7(0, t, x) > 0}. 

The analogous conjecture is that that this equation makes sense, that it has has a critical value L c , and that the 
limit *(/3)//3 2 exists and equals L c . 



7 Death 

7.1 Death in dimension d — 1 

We will show that death is certain when 7 = 0, /3 > and d = 1. 

Lemma 23. There are closed bounded intervals C (00, oo), indexed by fi 6 M.(R), and probabilities pp > 0, 
for j3 > 0, so that when /1 is compactly supported 

Q*f Q [u [Q , oo) {i{nY) = o)> P p. 

By arguing as in the proof of Lemma I12| one may check that 

Q^iu^in = o] = QlX'ffpM^) = 0] (7.1) 

whenever J is a closed interval supporting fj, and satisfying I C (— M, M). We define stopping times iteratively as 
follows. Let n = inf{£ : Ut(I(n) c ) > 0} (with inf{0} = +00). Supposing r„ < 00 we set 

r n+1 - inf{f > r„ : U t (I(U T J c ) > 0}. 

If t„ = 00 then we let t„ + i = 00. Set p(fi,f) — Q^'f'°[U[o.oc) (^(/ i ) c ) > 0]. The decomposition in Lemma[4]and 

(|7.1[) show that / — > p(fi, /) is non-decreasing. By the lemma above, Qai'°[ T i < 00 ] = 1) 5= 1 — By the 
strong Markov property at r n we have on the set {r n < 00} 

Q%l '°[rn+i < oo|W r „] = p({7 r „,exp(-{7(0,r„))) 
< p({7 r „,l) < l- P/3 . 

Then iterating the lemma shows that 

Q%i'°[U[a,oc) is not compactly supported] < Q^f '°[r n < 00] < (1 - pp) n 
completing the proof of death by the characterization in Lemma 1111 

Proof of Lemma 1231 Choose a closed bounded interval I = 7(/i) C (—00, 00) so that [^[o,oo)(^ c ) = 0] > 1/2. 
This is possible since paths die out and are compactly supported. Without loss we may assume I has length at least 
2. Now set I = 1+ [-1,1]. Choose M so that I C (-M,M). By $71$ we must show <3^ M,M) ' /3,1 [{/[ ,oo)(^ c ) = 
0] > P/3- Construct processes (u~ 7 u + 7 v + ) so that u~ is a DW((— M, M),0) process started at ^, and conditional 
on cr{u~} the process (u + ,v + ) is a solution to p.l[) on (— M, M) with parameters (/3, 0) and initial conditions 
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Uq = f3(l — e u (°'°°' x ^)dx and Vq — e u ( >°°>. Then Lemma [5] (i) shows that - + - has the same law 
as f/[o,oo) under q^ ] m > m )>' 3 > ( \ Q n ^ ne se ^ | u - ,(J C ) = 0} we have < /3I(x e I)da;. When this occurs, Lemma 
[6] shows that the total occupation measure it^ ^ becomes stochastically larger if we replace the initial conditions 
(u+,v+) by {(31{x € I)dx,I((-M,M)\I)). Therefore on { M[0 oo) (/ C ) = 0}, 

pi< oc) (i c (ri) = o\*{u-}} = Qi- + y w [u [0t00) (n = o] 



> 



u ,v 

(-M,M),P,Q 



f3I(xei)dx,I((-M,M)\iy ■ t ' 00 ^ I ~ J 

R t° e I W „i(M[%oo)(/ c ) = 0] byC3 



^ f3I(x£l)dx, 1(1") 

> exp(— /? Jw(x)dx) using (|3.17p 

provided that u> solves 

Aw < |w 2 — ?7W on int(J) and w(x) f oo as a; — > 9/ (7-2) 

with smooth 77 > satisfying r/ > f3 on I \ I. The proof will be complete if we find such a function u> and show 
that there is an upper bound on fjw(x)dx that is independent of the length of the interval /. 

Choose smooth non-increasing rj : [0, 00) — > [0, 1] so that rj(x) = 1 for x E [0, 1] and 77(33) = for x > 2. There 
is a unique non-increasing w(x) > 0, for a; € (0, 00) that solves 

Aw = |tti 2 — 2(3rjw on (0, 00), uJ(ir) — » 00 as x — ► and W(a;) — > as a; — * 00. 

Moreover the decay as x — > 00 ensures that L w(x)dx < 00. Suppose = [a, 6] and / = [a + 1,6 — 1], 

where 6 > a + 4. Then u>(a;) = w(a; — a) + w(& — a;) satisfies the required properties (17. 2|) with the choice 
77(3;) = /3T/(x — a) + ^77(6 — as). □ 

Remark. The extinction estimates for a DW(i?, 7) can be used to show that Q^ ) ' 7 [[/[o. 00 )([— L] c ) > 0] ~ 
exp(— 7 1 / 2 L) for large L. This in turn yields the finiteness of small positive exponential moments of J(l — 
exp(—2U(0,oo,x))dx which implies, via Novikov's criterion, that the exponential martingale t — > £ t (M^' 0,1 ) is 
uniformly integrable for small (3 and hence certain death. We omit the details since this simple method does not 
seem to be applicable for large (3 in d — 1 or any [3 > in dimension d > 1 . 



7.2 A block estimate for death 



By Lemmas [TT] and [12] to establish certain death for some parameter values (/?, 7) it is sufficient to choose a 
non-zero compactly supported \x and show that 

lim qfif" [E72?-=0]=l. (7.3) 

It can be simpler to verify the following block estimate, which will imply (|7.3p via a simple iteration reminiscent of 
the more direct branching process comparison used in [9]. 

Lemma 24. Suppose that there exists L,M > so that for all \x supported in Dl = [—L,L] d and satisfying 
/j(l) < M the following bounds hold: 

C A7 [U^ f 0] < 60 and Q°f>^ [U™»{\)] < e M, (7.4) 

where £q = 1/(4 ■ 3 d ). Then |7.ff[ ) holds for \i = MSq and hence certain death holds for 7). 
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Proof. We want to break a measure /i € M(R d ) into submeasures supported on blocks of size L and of mass at 
most M. We consider translates x + [-L, L) d for x <E 2LZ d and break the restriction of fi on x + [-L, L) d into at 
most Int[Af~ 1 /i(x + [-L, L) d )) parts, each of measure at most M (here Int[z] is the smallest integer greater than 
or equal to z). We require only that this splitting is done in a measurable manner. Let Nl^mIp) be the number of 
pieces into which /j, is thus decomposed, that is 

Nlm{^)= tok[M- l n(*+[-L>L) d )]. 

xeLZ d 

We now inductively define a sequence (X n : n > 0) of random measures. Let Xq = MSq. Let X\ = ii£} ,dD3L for a 
solution (u^\v^) to on D^i, with initial conditions = Xq and Vq = 1. The spatial Markov property 
(Lemma [7]) and the monotonicity of exit measures with respect to the initial nutrient level (Lemma [4] (ii)), imply 
that for R>3L 

Qxt;^ [uL Dr =o)>e [c%?f«[u™« = o] . 

Then inductively, provided that X n -\ =/= 0, we choose one of the submeasures X n _i that form the decomposition of 
X n _i as in the first paragraph (in some measurable manner), and denote its supporting cube as x n -% + [-L, L) d . 
Let D(n) = a;„_i + (— 3L,3L) d . Conditionally independent of (u^\v^ : j = 1, . . . , n — 1) we run a solution 
to (jl.ip on D(n) with initial conditions = X n _i and — 1. We then set 

X n = X n - l -X n ^ l +ut ) ' aD{n) . (7.5) 

When X n -i = we set X n = 0. 

We apply the decomposition Lemma [5] (i) to the domain D(n) to see that the total exit measure on dD{n) 
started from {X n _i\jji n ), 1) can be constructed in two parts as u^ ,dD ^ + u^ 1 ' 90 ^ where ii^ is a solution on 
D (n) started at -Xjj-ll-Dnv) ~ -^n-i and initial nutrient level exp(— u^ n \0, oo)). The monotonicity of exit measures 
in the nutrient level shows that u^' dDlj ^ is stochastically smaller than the total exit measure of a solution on 
D(n) started at (X n , 1). Noting that X n is certainly supported inside D(2n+i)L, the spatial Markov property then 
gives, for R > (2n + 1)L, 

n D R ,p n \tt8Dr 



'X A 



U%>* =0]>E Qj£f = 0] . (7.6) 



Using the simple inequality Int[z] < z + I(z > 0), the hypotheses of this lemma imply, when X n 0, that 

E[N LM {ut ] ' dD{n) )] < 3 d (e + e ) < 1/2. 

Note from (E3 that N LM {X n ) < N L ,M{ X n-\) - 1 + N LM {u ( ^' dD(n) ) when X„_ a f 0. Then n N LM {X n ) is 
a N valued supermartingale whose only possible limit point is zero, which implies that X n — for large n, almost 
surely. Then (|7.6p implies the sufficient condition (|7.3|) for certain death. □ 

As an illustration of this block condition, wc shall use it to show for any (3 > there exists 7 < j3 so that certain 
death occurs for (j3, 7). This implies that $(/3) < /3. The following lemma will imply that the block estimates will 
hold when 7 = /?, and a perturbation argument will ensure they hold 7 sufficiently close to (3. 

Lemma 25. Consider solutions in the domain D3 = (—3, 3) d to i/ie equation 

d t u = b-^-^^Au + b(3u(v - 1) + ^/uW, , 7) 

dtv = —uv, vq = 1. 

Then, given e > 0, £/iere exists bo = bo (/3, <i, e) > so £/ia£ /or b > bo, and whenever uo — f-i is supported in D\ and 
satisfies /i(l) < 1, 
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Remark. Note that (|7.7|) comes from (jl.ip via the scaling lemma with the choices e = 1, a = b, c = b 2 / d . Undoing 
this scaling shows that the block estimates (|7.4[) hold when 7 = j3 and L — b 2 ^ d , M = b with b = 6o(/3, d, eo)- 



Proof. We first consider solution to (|7.7[) on the smaller domain D2 and estimate the expected total exit measure. 
Using Lemma[6]we can and shall, at the expense of obtaining a stochastically larger exit measure, change the initial 
condition to /i + I(D 2 )(dx) and the initial nutrient to Vo = (1 — (6/3) _1 )/(£> 2 )- With these changes the reaction 
term b(3u(v— 1) is at most — u. Let <f>b. r (x) solve b^^~ d ^ d A(f>b tr = 4>b. r on the domain D r with boundary conditions 



(f)b. r = 1 on dD r . Calculus shows for r > 2 that u t ((f>b,r) + wf D2 {4>b,r) is a non-negative supermartingale. Taking 
expectations and r J. 2 we find £'[m^P 2 (1)] < J B (pb^^inidx) + dx). As & — > 00, (f>b^{x) decreases to zero as 
— * co for x G £>2 (argue, for example, using the probabilistic representation) and therefore uniformly for x € Di. 
This implies that E[u™ 2 (l)] -> as b — > 00, uniformly over /x as in the statement of the lemma. 

For a solution to (|7.7[) on Z?3 we apply the spatial Markov property Lemma [7] with the subdomain D^. Using 
also extinction estimates as in (I3.17P this shows that, when fj, supported in D\, 



E[u^(l)]<E[u^(l)] and P[u<£° (1) > 0] < E[l - 



provided that b^^^^Aw < w 2 /2 on D 3 and M{w{x) : x e dD r } — > 00 as r ] 3. The test function u>(a;) = 
12^=1^ + 3)~ 2 + (3 — Xi)~ 2 satisfies these requirements provided b > 1 and also sup{w(a;) : a; S SZ^} < 00. 
The lemma now follows from our control on the expected exit mass -E[u^P 2 (l)]. □ 

The proof that \&(/3) < /3 will follow from Lemma[M]once we have shown that on a box D, the block estimates 
Qft'i ' 1 [^oa ) 7^ 0] anc * Q^i^lE^PCl)] are continuous as 7 f /3, uniformly over /x supported in a certain strict sub- 
box and with a certain bounded total mass. This follows by a change of measure argument. Indeed the derivative 



dQ^ 1 y /dQn 1 on i s given by the exponential martingale £ t ((/3 — j)M(l)) where M t (l) is the martingale part 
of U t (l). Note that [M(l)] t = f7 [0)t] (1). The uniform integrability of £ t ((f3 - j)M (1)), when (3 - 7 is small, follows 
from the finiteness of the exponential moments 



[exp(AC/ [0 , 00 )(l))] 



< 



yD.o 



[exp(AE/ [0)Oo) (l))] =e^ ( 



< OO 



where <f>^ solves, for small enough A > 0, = -(1/2)(0( A )) 2 - A on D and cj> { ^ = on clD. Then 



Q?,f 7 [t^(i)]-02f 7 [uL D (i)} 



< (o£f [o* B (i)(«»W-7W)) - i)f 

< [(^ D (l)) 2 ] 0?f ^ [(foo((/3-7)M(l)) - 1)= 



Then Q°f' [(U^ D {l)) 2 } < Q°'°[(U^ D (1)) 2 } can be bounded in terms of the total mass while, when - 7 is 
small enough, 



< 



[(£«,((£ -7)M(l))-l) a ] 

e 2(/3- 7 )M oo (l)-(/3- 7 ) 2 (7 [0iOo) (l) 
1/2 



= (cxp(^ (A) )) 



[f 00 (4(/3- 7 )M(l))] 
1/2 



,6(/3-7) 2 !7 [0l «o(l) 



1/2 



e 6(/3-7) 2 (7 [0l oo)(l) 
1/2 



1 



- 1 



where A = 6(/3 — j) 2 . Since decreases to as A J, 0, this shows the required continuity in 7. The argument for 
Q^P™ * 0] is similar. 
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7.3 Proof that *(/3) < c x [3 2 in d = 3 



We sketch a rather simple argument for death, based on the decomposition in Lemma O (i). This decomposition 
suggests we may construct the total occupation measure for a solution to (jl.ip in two parts: * from a pro- 

cess u~ with law ' 7 , and » from a solution u + to that conditional on a{u~} has initial condition 

Uq(cIx) = /3(1 — exp(— u~(0, co, x))<ia;. The expected initial mass can be exactly calculated via the Laplace 

functional of u~(0, oo, x) and shown to be bounded by co/37 _1,/2 £t(l). This suggests, when co/?7 -1 / 2 < 1, that this 
decomposition can be iterated leading to a convergent geometric series and a finite total occupation measure. 

We now formalise this argument, working on finite domains and with exit measures. 
Lemma 26. When c /?7~ 1/2 < 1 we have 

oo 

Q D ^ [c£ D *(i)] < M0 (L) ) + M(i)E( c ^7' 1/2 ) fc 

k—1 

where <p solves Acp^ — j(f>^ on Dl and <jy^ = 1 on ODl, and Co = J 3 ijj{x)dx where (ip(x) : i? 3 \ {0}) is the 
maximal solution to 

Atp = \%p 2 + ip on R 3 \ {0}, %/}(x) -> as \x\ -> oo and V(^) -> oo as |x| -> 0. (7.8) 



Remark. The maximal solution ip can be constructed as follows: non-negative solutions ip( 6 ' to Aip^ = 
(l/2)(ip^) 2 -\-tp( e > on {x : < e < |x| < oo}, with boundary conditions ip^\x) — ► as |x| — ► oo and ip( e '(x) — * oo as 
|a;| | e exist and are unique (take limits through increasing boundary conditions, using the upper bound C(|a;| — e)~ 2 , 
for existence and use the maximum principle for uniqueness.) Moreover ip^ decrease as e J, to function ip de- 
fined on R 3 \ {0}. Interior regularity estimates show that that ip solves (|7.8p . The boundary condition for ip^ at 
|a;| = e imply, by a comparison argument, that ip is the maximal solution. Finally, another comparison shows that 
< ip < 4| 2 and that ip has exponential decay at infinity so that Cq = J ip(x)dx < oo. 

Proof. L,(3,j are fixed throughout. For convenience in this proof we let F(fi,f) — Q I ^ L f^ n \U^ )L (1)] and let 
S(/z,aV) be the the measurable kernel given by the law of [3(1 — exp(—U(Q,oo,x))dx on A4(Dl) under Q^ Ln ■ 
Then the decomposition in Lemma [5] (i) , monotonicity and the exact formula for exit measures of superprocesses 
(derived, say, from (|3.8p ) imply 

F(ji, 1) = Q^i [ug>r. (i)] + Q^n _ e -u(o,oo,x) dXi e -tr(o,oo,*)y 

< Q° Ln [u™ L (i)] + Q° Ln \f(P(i - e- u ^°°^dx, 1) 



(7.9) 



Define the range by 72 — Us>o Ut>s support(LTf), where support(/i) is the closed support of fi. Then 



l_ e -^(o, 



dx 



< Qu 1 " 7 [W e 72.] = /3 / n(^' x ^)dx 



where ?/> = t/^ 7 ' 21 ) defined for y £ Dl \ {x} is the maximal solution to 

Atp = ^ip 2 + jip on Dl \ {x}, ip = on ODl and ip(y) — > oo as |y — x| — > 0. 

The final equality follows using mostly the arguments in [10] Theorem III. 5. 7 (a) (and constructing tp^'^ as in 
the above remark). A comparison argument shows that ip^ 1,x \y) < r yip(^f 1 ^ 2 (y — x)) and this leads to the bound 
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/i/(l)H(/i,di/) < co^-^VCl)- Note that </>( L ) < 1. We can now interate the inequality (|7.9p and the lemma 
will follow once we have show that the remainder term converges to zero, and this will follow from the fact that 
F(n, 1) — > when — » 0. To see this, note that 

FGM) < Q^[Ui D Hl)I(r>t)]+Q^[u^(l) 

...1/2 



< 



{Q D ^[{UL D ^l)) 2 ]Q D ^[r>t]) +Q^[urni) 
Sfe)* [(^(l)) 2 ] ( ?m,i' /3 ' 7 [* > *]) V2 + 



where 86 = (/3 - j)6 on [0,t] x D L , O = and 6 = 1 on [0,i] x <9D L . Note the term Q°+p I{DL)dx [{U^ L {1)) 2 } 

is bounded by C(L,(3,j) < oo for /i(l) < 1. The term Q^\ ,t3,1 [r > t] converges to zero as £ — > oo, uniformly 
over < 1, for instance by the estimate in Lemma UM These together allow us to see that F(p, 1) — * when 
ju(l) — * completing the proof. □ 



To complete the proof that \&(/?) < ci/3 2 we will apply the lemma above to a sequence of domains D(n) 



D, 



i , and show that the expected exit measures decay geometrically. Scaling shows that the solution 



(7- 1/2 ) 



from the above lemma satisfies 1 \x) — ifi^^x) where A(f> = cj) on D\ and = 1 on 8D\. Moroever a 
comparison argument shows that 1 2 \x) < </>( 7 1 2 -*(0) = </>(0) < 1 for x S 8D(n). The lemma therefore 

implies that 



(7.10) 



whenever /i is supported on 8D{n). We may now choose 7 = ci/3 2 with ci < 00 large enough that the right hand 
side of (|7.10p is at most (1/2)(1 + -0(O))/i(l) < m(1). Iterating shows when fj, is supported on dD(l) that 



n D(n+l),/3,c 2 /3 2 



c/ ^(n+D ( i ) l <((i + 0( O ))/2yv(i). 



This implies certain death. Indeed the spatial Markov property and the extinction estimate (I3.17|) show that 



where 



U°?«> ? 
)=£(V 



< 



1 - e 



as n — > 00 



12 

(xj + n7~ 1 / 2 ) 2 



2/3- 



12 



(n.7 - 1 / 2 — Xi) ; 



satisfies Ait; < (l/2)w; 2 — (3w on D(n). Certain death follows from Lemmas fTTI and [T2l 



Remark The approach above can be applied to the case of d = 2 and small (3. An analysis of the first moment 
shows that when 7 = exp(— C//3) the above series construction converges for suitable C and hence certain death. 
We do not include the details since the argument is quite crude and it seems easier to conjecture that there is death 
when 7 = for small (3. Note however that when 7 = a first moment argument will not show death since the 
first moment Q?i '°[C^P R (1)] — * 00 as R — * 00. We hope to comment on this in a subsequent paper. 
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